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O . ABSTRACT 
(N 

Context. Studies of the generation and assembly of stellar populations in galaxies largely benefit from far-IR observations, consider- 
ed . ing that the IR flux is a close prior to the rate of star formation (the bulk of which happens in dust-obscured environments). At the 
same time, major episodes of nuclear AGN accretion are also dust-obscured and visible in the IR. 

Aims. At the end of the Spitzer cryogenic mission and the onset of the Herschel era, we review our current knowledge of galaxy 
evolution at IR wavelengths, and model it to achieve as far as a complete view of the evolution of cosmic sources. We also develop 
new tools for the analysis of background fluctuations to constrain source counts in regimes of high confusion, as it happens for the 
Herschel sub-mm surveys. 

Methods. We analysed a wide variety of new data on galaxy evolution and high-redshift source populations from Spitzer cosmologi- 
cal surveys, and confront them with complementary data from mm ground-based observations and constraints from the far-IR diffuse 
radiation, as well as preliminary results from Herschel surveys. 

Results. These data confirm earlier indications about a very rapid increase in galaxy volume emissivity with redshift up to z ^ 1 
(— I ' [p(z) oc (1 -I- z)*], the fastest evolution rate observed for galaxies at any wavelengths. The observed Spitzer counts require a combina- 

tion of fast evolution for the dominant population and a bumpy spectrum with substantial PAH emission at z ~ 1 to 2. Number counts 
I . at long wavelengths (70 through 1100 pm) confirm these results. All the present data require that the fast observed evolution from 

O ' z=0 to I flattens around redshift 1 and then keeps approximately constant up to z ^ 2.5 at least. Our estimated redshift-dependent 

^ ' bolometric comoving energy density keeps lower at z ^ 1.5 than some previously published results based on either large extinction 

C/j I corrections, or large spectral extrapolations. 

, Conclusions. The present-day IR/sub-mm data provide evidence of a distinct population of very luminous galaxies becoming dom- 

inant at z > 1. Their cosmological evolution, peaking around z - 2, shows a faster decay with cosmic time than lower luminosity 
systems, whose maximal activity is set around z - 1, then supporting an earlier phase of formation for the most luminous and massive 
galaxies. From a comparison of our results on the comoving IR emissivity with recent estimates of the redshift-dependent stellar mass 
functions of galaxies, we find that the two agree well with each other if we assume standard recipes for star formation (a universal 
Salpeter IMF) and standard fractions (~ 20 - 30%) for the contribution of obscured AGN accretion. Systematic exploitation of the 
forthcoming Herschel survey data will be important for confirming all this. 
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0^ 1. Introduction activity (Kormendy and Sanders 1992; Franceschini et al. 1994; 

' Sanders and Mirabel 1996). Mid- and far-infrared observations 

Galaxy evolution has been studied in detail and is currently quite of the high-redshift universe are then needed for precise quantifi- 

. _ precisely known between redshift and approximately 1 , mostly nation of the dust-extinction effect: without a direct observation 

^ based on deep optical imaging and spectroscopy (see, among of the dust-reradiated pait of the young stellar population emis- 

5^ others projects, DEEP-2, Faber et al. 2007; COMBO-17, Bell ^ion, a reUable estimate of the start formation rate is difficult, 

. C3 , et al. 2004; K20, Cimatti et al. 2002; COSMOS, Scoville et al. or even impossible (e.g. Poggianti & Wu, 2000; Poggianti et al. 

2007; VVDS, Le Fevre et al. 2005). One important aspect of 2001) 
these investigations concerns the fraction of optical-UV radiant 
energy emitted by young stellar populations, which is absorbed 

by dust in the ISM medium and re-emitted at long wavelengths. the redshift interval from local up to z ~ 1, relatively 
It has been shown that, while being modest for evolved galax- deep IR observations have been obtained with early space in- 
ies like early-type spirals, and negligible in E/SO galaxies, this fr^red missions like IRAS (Soifer et al. 1988) and ISO (Genzel 
fractional energy may become strong or even dominant during & Cesarsky 2000). Based on the fair sensitivity and moderate 
active phases of star formation, galaxy assembly, and nuclear spectral coverage of the latter, the galaxy and AGN evolution in 
the redshift interval < z < 1 has been mapped (Franceschini 
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For a long time, both theoretical expectations (Franceschini 
et al. 1994, 1997; Blain & Longair 1993; see also Kormendy 
& Sanders 1992) and preUminary observational results by mil- 
limetric and sub-millimetric surveys with SCUBA (Smail et al. 
2002; Hughes et al. 1998; Blain et al. 1999, 2002; Chapman et 
al. 2005) and IRAM (Greve et al. 2005; Bertoldi et al. 2007) 
have produced evidence that major events of star formation in 
the universe happened at high redshifts, z > 1, and are likely re- 
sponsible for the origin or assembly of the most massive galaxies 
today. Such events would be accompanied by major energy re- 
lease, but, because of heavy dust extinction, would be detectable 
only through far-infrared observations. Sub-miUimetric surveys 
indeed have identified sources that are extremely luminous at 
long wavelength, but almost undetectable in the optical (e.g. 
Chapman et al. 2005; Mobasher et al. 2006; Rodighiero et al. 
2007). 

A dramatic improvement in our knowledge of the infrared 
emissivity of source populations at redshifts z > 1 , correspond- 
ing to the the era of peak star formation and AGN activity, 
has recently been made possible with the deep sky exploration 
by the Spitzer Space Telescope (SST, Werner et al. 2004). The 
Multiband Imaging Photometer (MIPS, Rieke et al. 2004) on 
SST combined sensitivity and spectral coverage, particularly 
with the Aej f = 24 jim band, adequate for identifying and char- 
acterizing substantial numbers of far-IR galaxies up to z ^ 3 (Le 
Floch et al. 2005; Dole et al. 2006; Caputi et al. 2007; Lonsdale 
et al. 2004; Rowan-Robinson et al. 2005). The MIPS instrument 
also includes photometric imagers at A^fj - 70 and 160 jjm (al- 
though the latter was limited by source confusion). A further cru- 
cial contribution by the SST mission was the capability ofi'ered 
by the Infrared Spectrometer (IRS, Houck et al. 2004) to obtain 
high-quaUty far-IR spectra of high-redshift sources for the first 
time. 

A historical difficulty in identifying and interpreting lumi- 
nous high-redshift sources, i.e. the diflicult association of the op- 
tical counterparts to the IR-mm sources due to the poor spatial 
resolution and the large IR error-boxes, is also largely overcome 
by observations with the Spitzer Infrared Array Camera (IRAC, 
Fazio et al. 2004). IRAC deep multiwavelength imaging at 3 to 
10 /im allows continuous band-merging and association of far- 
IR and optical/near-lR source catalogues in this critical spectral 
transition region (e.g. Egami et al., 2004). 

The multiwavelength observations by Spitzer have fully con- 
firmed previous expectations of numerous populations of lumi- 
nous galaxies and active nuclei at high redshifts. Sources at z > 1 
and up to ~ 3 are routinely detected with the sensitive, high- 
spatial resolution 24 /im channel of Spitzer/MlPS (Lacey et al. 
2004; Perez-Gonzales et al. 2005; Lacey et al. 2004; Caputi et 
al. 2006; Daddi et al. 2007). The MIPS imager benefits from a 
wide enough detector format to allow deep surveying with high 
mapping speed of several areas of relevant cosmological inter- 
est (Papovich et al., 2004; Lonsdale et al. 2004; Le Floch et al., 
2004; Shupe et al. 2007; Sanders et al. 2007). 

MIPS surveys at longer wavelengths (70 and 160 yum) have 
produced more limited samples of distant sources (Frayer et al. 
2006a,b; 2009), because of the much longer integration times 
required and the lower mapping speed at faint flux levels. These 
faint far-IR Spitzer sources are also more difficult to identify, 
because of the poor Spitzer resolution at such long wavelengths. 
A major step forward is now possible with the Herschel Space 
Observatory (Pilbratt 2005; 2010). 

The present paper is dedicated to a statistical analysis 
of a very extensive dataset of IR sources, aimed at deeper 
understanding of the cosmological significance of the long- 



wavelength emission in the complex process of galaxy formation 
and evolution. A previous analysis by Franceschini et al. (2001, 
AF2001), based mostly on IRAS and ISO data, as well as on es- 
timates of the cosmic IR background (CIRB) by COBE, has al- 
ready achieved interesting results. These concerned in particular 
the evidence of a population of fast-evolving starburst galaxies 
between z=0 and 1 and the indication of a constant IR emissivity 
for galaxies at z > 1 . Results consistent with these were reported 
by Chary & Elbaz (2001) and Elbaz et al. (2002). Another in- 
dication by AF2001 was in favour of a top-heavy stellar initial 
mass function (IMF) during major star-formation phases, to rec- 
oncile the high IR emissivity with local contraints on the galaxy 
stellar mass function. 

All these had to be taken as preliminary hints, owing to the 
lack of any statistical data above z ~ 1 available to AF2001. A 
much more detailed and complete analysis is now possible based 
on the vast amount of new information available from Spitzer, 
particularly at high redshifts. New essential data also include ex- 
tensive sub-millimetric imaging surveys with BLAST (Pascale 
et al. 2008), and preUminary constraints by the Herschel obser- 
vatory. Thanks to the tight relationship of far-IR emissivity of 
galaxies with the instantaneous rate of star formation (Kennicut 
1998; Poggianti et al. 2001 ; Elbaz et al. 2002), a model reproduc- 
ing these observations of faint IR galaxies will provide us with a 
detailed description of the history of star formation in galaxies. 
While other galaxy evolutionary models have recently been pub- 
lished (Lagache et al. 2004 and references therein), the current 
paper accounts for a lot of new, partly unpublished information, 
particularly on the redshift and luminosity distributions of faint 
Spitzer/MIPS sources, time-dependent luminosity functions, and 
background fluctuations, providing us with an accurate descrip- 
tion of the high-redshift evolution. 

The present paper also includes a first quantitative compari- 
son of the history of the bolometric IR emissivity with the evo- 
lutionary stellar mass functions of galaxies (from deep Spitzer 
IRAC near-lR surveys), entailing interesting constraints on the 
nature of the high-redshift luminous source population. We re- 
call that, while the CIRB background intensity is a record of the 
past star-formation activity biased to low redshifts by the usual 
redshift factor (1 + z)"' (e.g. Harwit 1999), the comoving mass 
density in galactic stars is a completely unbiased tracer that is 
much more sensitive to the high-redshift activity. 

The paper is structured as follows. Section 2 reports on the 
wide variety of observational data on distant IR sources mak- 
ing our reference dataset. These include some new estimates of 
source number counts from the SWIRE surveys at 70 and 160 
fitn and new evaluations of the redshift distributions of faint 
24 pm sources. We also report new data here on local lumi- 
nosity functions in Spitzer/MIPS bands, establishing a bench- 
mark against which evolutionary models of cosmic sources are 
compared. In Sect. 3 we report on multi-wavelength models of 
galaxy and AGN evolution fitting the whole variety of exist- 
ing data, including the local background spectral intensity, as 
well as data on distant and high redhift sources. In Sect. 4 we 
discuss our results and identify particular evolutionary patterns 
emerging from our analysis, including detailed comparison of 
the source comoving IR emissivity with the comoving stellar 
mass density and discuss implications from this comparison. Our 
conclusions are summarized in Sect. 5. Two appendices are ded- 
icated to discussing the contribution of discrete sources to the 
background intensity and cell-to-cell fluctuations and to elabo- 
rating some model predictions, with particular reference to the 
forthcoming Herschel observations and source confusion linaita- 
tions that will be encountered. 
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Fig. 1. Euclidean-normalized differential number counts of ex- 
tragalactic sources at 24yum compared with our model fit. The 
red circles are from the analysis of SWIRE survey data by 
Shupe et al. (2008), black squares from Papovich et al. (2004). 
The contribution by type-I AGNs is shown as green dot-dashed 
line, moderate-luminosity starbursts (the LIRGs) make the cyan 
short-dash line (type-II AGNs and starbursts are included in the 
same population on the assumption that in both classes the IR 
spectrum is dominated by starburst emission). The red long- 
dashed line corresponds to the population of high-luminosity 
sources dominating the IR emissivity at high redshifts. The dot- 
ted line is the separate contribution of normal spirals, while the 
continuous line is the total model counts. 



All quantities are computed in this paper assuming a geom- 
etry for the Universe with //q - 70 km s"' Mpc"', Q,„ - 0.3, 
Qa = 0.7. We indicate with the symbols 5 24 and L24 the flux 
density in Jy at 24 fj.m and the monochromatic luminosity vLy at 
24 yum in erg/s normalized to the solar value Lq (and similarly 
for other wavelengths). 



2. The dataset 

Our analysis reported in the present paper has benefited from 
the availability of highly reliable deep extragalactic surveys in 
a variety of IR bands, thanks mostly to the sensitivity and pho- 
tometric accuracy of the Spitzer Space Telescope's MIPS instru- 
ment (Rieke et al. 2004). These properties of the far-IR cam- 
era have allowed not only for probing to very deep flux limits 
(S 24 - 20 fiJy) in the 24 /im channel, but also complementing 
these ultradeep maps with large-area surveys covering several 
tens of square degrees, still to appreciable depths. 

The main Spitzer's limitation rests in the poor spatial reso- 
lution offered by the diffraction-limited photometric imagers in 
the critical far-IR (70 to 160 //m) spectral region, which implies a 
strong confusion noise, particularly at 160 jjm. This limit is cur- 
rently being overcome by the ongoing operation of the Herschel 
Space Observatory. We are also recovering important constraints 
on the high-redshift galaxy emissivity at long wavelengths from 
deep millimetric and sub-mm observations with large ground- 
based observatories (JCMT, IRAM, CSO among others). 



Fig. 2. Euclidean-normalized differential number counts of ex- 
tragalactic sources at 70/im compared with our model fit. Red 
open circles are from our analysis of the SWIRE survey data 
(Vaccari et al. in preparation, see text), open squares to a very 
deep survey by Frayer et al. (2006a), the cyan circles are the most 
recent assessment of the counts based on Spitzer/MIPS observa- 
tions in the COSMOS area by Frayer et al. (2009). The red filled 
datapoint at bright fluxes is from the IRAS 60 jum survey appro- 
priately scaled to the current effective wavelength (Vaccari et al. 
2009). Moderate-luminosity starbursts (including type-II AGNs) 
make the cyan short-dash line. The red long-dashed line corre- 
sponds to the population of high-luminosity sources dominating 
the IR emissivity at high redshifts. The dotted line is the sepa- 
rate contribution of normal spirals. At bright 70 yum flux densities 
type-I AGNs, shown as green dot-dashed line, provide an impor- 
tant contribution to the counts. The continuous line is the total 
model counts. 
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Fig. 3. Euclidean-normalized differential number counts of ex- 
tragalactic sources at 160jum compared with our model fit. Lines 
are as in Fig. |2] The red open square datapoints are from a re- 
analysis of flie Spitzer/MIPS data by Bethermin et al. (2010). 



4 



Franceschini, A., et al.: Galaxy Evolution at IR Wavelengths 



2.1. The IR multi-wavelength source number counts 

2.1 .1 . The Spitzer MIPS 24 /um counts 

The Spitzer Wide-Area Infrared Extragalactic (SWIRE) survey 
(Lonsdale et al. 2003, 2004), the largest of the Spitzer extragalac- 
tic Legacy Science programmes, mapped 49 square degrees split 
into 6 widely separated fields, at wavelengths from 3.6 to 160 
microns using both IRAC and MIPS. The MIPS 24 jjm counts 
were derived by Shupe et al. (2008) from all six SWIRE's areas 
to minimize and averaged out sample variance effects. Indeed, at 
fluxes brighter than 5 24 - 2 mJy, large field-to-field variations, 
by up to a factor of two, are found in the counts, whereas they 
showed differences on the order of only 10% in the sub-mJy re- 
gion. The Euclidean-normalized differential counts reported by 
Shupe et al. (2008) averaged over the six independent SWIRE 
fields are reported in Figure [1] 

Deeper extragalatic number counts, down to flux limits of 
5 24 - 20 fxJy with the MIPS 24 yum channel, have been obtained 
in the GOODS, GTO and COSMOS surveys by Papovich et al. 
(2004), Chary et al. (2004), Marleau et al. (2004), Rodighiero 
et al. (2006), Sanders et al. (2007). A selection of these deeper 
counts is reported in Fig. [T] It is interesting to note the perfect 
agreement among the results of different teams, sky areas, and 
data reduction procedures required corrections to the total fluxes 
of only few % altogether, showing the high reliability of MIPS 
at these wavelengths. 

Altogether, the number counts show a steady super- 
Euclidean slope as a function of flux density all the way down 
to the SWIRE completeness limit of 0.4 mJy, with a hint of con- 
vergence below it. At fainter fluxes the counts from the deeper 
surveys show a very fast convergence down to the MIPS confu- 
sion limit, estimated to fall around a few tens of fiJy. A similar 
bumpy behaviour, though not as much extreme, was observed by 
Elbaz et al. (1999) in the deep ISO counts at 15 pm and was 
interpreted as a signature of very fast evolution in the galaxy 
emissivity from z=0 to 1 (see Sect. 3). 

2.1 .2. The Spitzer MIPS counts at 70 and 1 60 pm 

MIPS observations at longer wavelengths required much more 
extensive integrations to achieve comparable depths. Deep sur- 
veys have then only been possible on small sky areas, the bulk of 
the observations in the SWIRE and ELS fields having been per- 
formed to relatively shallow depths. The deepest observations at 
70 and 160 pm were also limited by source confusion. 

We used the SWIRE database to estimate the extragalactic 
number counts at 70 pm. While a detailed account of this work 
will be reported in a separate paper (Vaccari et al., 2009, and pa- 
per in preparation), we summarize here the main facts. Source 
catalogues at 24 and 70 pm were extracted by us in the six areas 
observed by the SWIRE survey over a total of 49.5 square de- 
grees, where good complementary data exist to perform an accu- 
rate subtraction of stars. As a consequence of the quite different 
sensitivities, the vast majority of the 70 pm sources had counter- 
parts at 24 pm. In turn, essentially all the 24 pm objects are iden- 
tified in the short wavelength IRAC maps. This has prompted 
us to obtain highly reliable catalogues at these wavelengths. In 
all, Vaccari et al. identified about 10000 extragalactic sources at 
70 pm brighter than the 5cr limit, respectively. The survey com- 
pleteness and source reliability were extensively tested by in- 
serting simulated sources into the real maps. The corresponding 
Euclidean-normalized differential counts from the SWIRE sur- 



vey are reported in Figs.|2] These data compare well with those 
derived by Frayer et al. (2006b) from the ELS survey. 

Deeper number counts at 70 pm, based on a very long MIPS 
integration in the central lO'xlO' area of GOODS-N and a shal- 
lower GTO observation over 0.6 sq.deg. in GOODS-N, are re- 
ported in Frayer et al. (2006a) and shown in Fig. |2] A survey 
at 70 and 160 pm on a much wider COSMOS 2 square de- 
gree area is reported by Frayer et al. (2009). These observa- 
tions have prompted them to accuratele evaluate the extragalactic 
source confusion noise at 70 pm, which is estimated to be around 
0.35 +0.15 mJy, and correspondingly estimated the counts down 
to a Act flux limit of 1 .2 mJy. 

We report in Fig.[3]the results of a systematic re-analysis of 
Spitzer/MIPS number counts data by Bethermin et al. (2010). At 
160 pm, the Spitzer Icr confusion noise sets around cr con fusion ~ 
9.4 + 3.3 mJy according to Frayer et al. (2009) and Dole et al. 
(2006, see also Sect.|4]below). The flux limit of ~ 40 mJy in the 
number count data of Fig.[3]is then close to that set by confusion. 
This issue of the source confusion noise for the Spitzer/MIPS 
observations will be discussed later in Sect. l2.2l when talking of 
background fluctuations. 

Interesting features can be found in Figs. |2] and [3] The 70 pm 
ultradeep counts reveal fast convergence at faint fluxes, although 
with limited statistics, similar to those observed at 24 pm (Fig. 
[U and 15 pm (AF2001). This implies that the majority of the 
extragalactic background light (EBL) has already been resolved 
into sources by Spitzer at this wavelength and that a major part 
of the radiant energy produced by sources has been identified. In 
particular, about 60% of the total estimated EBL at 70 pm has 
been resolved into sources by Frayer et al. (2006a) down to the 
confusion limit, and as much as 75% of it at 24 pm. 

Also the slopes of the differential counts in 'the super- 
Euclidean' part at moderately bright fluxes up to the peak bring 
interesting information. These counts are flatter at 70 pm and 
steeper at both 24 and 160 pm, as a consequence of the very dif- 
ferent K-corrections for the average galaxy spectra at the respec- 
tive wavelengths: the K-correction is minimal at 70 pm owing to 
the very steep spectra towards short wavelengths. This also en- 
tails different average source redshifts for the flux-limited sam- 
ples at the two wavelengths. 

2.1 .3. The CSO, SCUBA, and IRAM counts at 350, 850, and 
1200yum 

Since the early observations with SCUBA (Ivison et al. 1998), 
later confirmed with IRAM and CSO, deep surveys with milli- 
metric and sub-millimetric observatories have contributed some 
fundamental insight into the high-redshift universe, by show- 
ing the existence of extremely luminous source populations that 
were essentially undetected in the optical. The identification 
of the optical counterparts of these objects has turned out to 
be a painful process, because of the large positional error-box 
of mm telescopes and the source faintness (Blain et al. 2002). 
Eventually, such a process required a simultaneous identifica- 
tion in the radio (Chapman et al. 2005), or the high-resolution 
imaging with mm interferometers (e.g. IRAM-PdB, Tacconi et 
al. 2006), which are available for few sources. The situation for 
astronomy at millimetric wavelengths will not change signifi- 
cantly until the advent of ALMA. 

The limited information about statistical properties, like the 
high-z luminosity functions or redshift distributions of the mm 
sources, led us to mostly use the information in the number 
counts to constrain our model at long wavelengths (hence high 
redshifts, due to large K-correction). Figure |4] compares the re- 
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Fig. 4. Left: a collection of observed number counts at 350 jum by Khan et al. (2007) compared with our best-fit model prediction. 
Center: the observed number counts at 850 //m by Blain et al. (1999) and Barger, Cowie and Sanders (1999, open squares). Filled 
squares are a recent assessment by Coppin et al. (2006). Right: number counts at 1 100 /vm by Austermann et al. (2009) from an 
Aztec survey at JCMT in the COSMOS area. Line types as in Fig. |2] 



suits of our model with the integral counts at 350, 850, and 1 100 
jum. A collection of the 350 //m counts is reported by Kahn et al. 
(2007) based on deep CSO observations. 

Number counts at 850 jt/m with SCUBA on JCMT (central 
panel of Fig. |4| have been reported by Blain et al. (1999) and 
Barger, Cowie, and Sanders (1999). The most statistically accu- 
rate recent evaluation has been published by Coppin et al. (2006) 
based on the largest, most complete and unbiased extragalactic 
submillimetre survey, the SCUBA HAlf Degree Extragalactic 
Survey (SHADES). The actual survey covered 720 arcmin^ to 
2 mJy Icr, detecting approximately a hundred sources. The dif- 
ferential number counts by Coppin et al. are reported in Fig. [13] 
below. 

Counts at millimetric wavelengths (1100-1200 ;um) from 
IRAM MAMBO and BOLOCAM have been published by vari- 
ous authors. The right panel of Fig. |4]reports a recent assessment 
of the number counts at 1100 /im by Austermann et al. (2009) 
from a survey with Aztec at JCMT in the COSMOS area. 



2.1.4. The BLAST sub-millimetric observations 

The balloon-borne BLAST mission (Pascale et al. 2008) has re- 
cently performed an exploratory survey of a region in CDFS us- 
ing a similar instrumental setup and the same waveband filters 
that are used by the Herschel Space Observatory at these wave- 
lengths. Devlin et al. (2009) and Patanchon et al. (2009) estimate 
source number counts based on a P(D) analysis from deep sur- 
vey data in CDFS (0.8 sq.deg.) combined with shallower data 
over 8.7 sq.deg. A comparison of the BLAST estimated counts 
at 250, 350, and 500 //m with predictions of our reference model 
(Sect. [3]below) is reported in Fig. |5] 



The 

imaging 
imposinj 
wavelen: 
ground 
consider 
BLAST 
counts. 



BLAST survey data provide us with the first extensive 
survey of the sky at sub-millimetric wavelengths, hence 
I an essential constraint on the source emissivities at 
gths close to the peak of the CIRB extragalactic back- 
The uncertainties related to source confusion led us to 
the faintest datapoints at 10 to 20 mJy at the three 
wavelengths as upper limits on the real extragalactic 



TI, I 



S50 fj^m 



350 ^zm 




500 /zm 



100 
S [mJy] 



Fig. 5. Data on the Euclidean-normalized differential number 
counts of sources at three sub-millimetric wavelengths obtained 
by Patanchon et al. (2009, see also Devlin et al. 2009) from ob- 
servations with BLAST in the CDFS area. The two BLAST 250 
jim datapoints at the faintest fluxes have been combined in our 
plot. The black continuous line is the total predicted counts by 
our reference model. 



2.1.5. Preliminary constraints from Herscliel observations 

After its successful launch in May 2009, instrument validation 
and Performance Verification, the Herschel Space Observatory 
carried out Science Demonstration observations during the 
second semester of 2009 and has subsequently started the 
routine observation phase. Herschel carries on-board two 
sensitive photometric imaging cameras, PACS and SPIRE, 
covering a wide wavelength interval between 70 and 500 
jim, where the bulk of the CIRB energy resides. While a 
first systematic account of cosmological observations with 
Herschel will be reported in a dedicated issue of this Journal, 
some constraints can already be inferred from the Science 
Demonstration Phase Workshop in December 2009 (see 
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http : //her schel . esac . esa . int/SDP JR.wkshop . shtml). 

At the present stage we have confined ourselves to checking 
that the data appearing in those reports (particularly those by S. 
Oliver and D.Lutz) do not conflict with our model results. At the 
same time, these data appear to agree closely with the BLAST 
results at the corresponding wavelengths. 




Deflection [MJy/ sr] 

Fig. 6. Histogram of the distribution of background intensity in 
MJy/sr measured inside circular apertures of 12, 14, and 16 arc- 
sec radius from the FIDEL 70 jjm map (Frayer et al. 2009, green, 
cyan, red histograms respectively). The observed histograms of 
deflections per unit sky area Dg are compared with a Gaussian 
distribution representing the instrumental noise, having an rms 
of 0.62 mJy/beam (blue line), and with the probability distribu- 
tion of deflections Da due to sources, which are modelled ac- 
cording to our best-fit source-count model (red line). The black 
line is a convolution of the two and provides a good fit to the 
observational distribution. 




Fig. 7. Data on the cosmic infrared background (CIRB; e.g. 
Hauser et al. 1998) spectral intensity compared with the best-fit 
model predictions. The three lower datapoints in the far-IR are 
from a re-analysis of the DIRBE data from the all-sky CORE 
maps by Lagache et al. (1999), the shaded areas from Fixsen et 
al. (1998, black shade) and Lagache et al. (2004, red shade). The 
two mid-IR datapoints at 15 and 24 fim are the resolved fraction 
of the CIRB by Elbaz et al. (2002) and Papovich et al. (2004). 
The triangle datapoint at 8 jum comes from the integration of the 
IRAC galaxy counts by Fazio et al. (2004) and Franceschini et 
al. (2008), while the green triangle at 70 fim is from the ultra- 
deep Spitzer survey by Frayer et al. (2006) (recalibrated by a 
scaling factor of 1.2 as indicated by Frayer et al. 2009). The dot- 
ted histograms are upper limits on the CIRB intensity set by Te V 
cosmic opacity measurements (Stanev and Franceschini 1998). 
The thick line is the predicted CIRB spectrum from the multi- 
wavelength reference model for galaxy evolution discussed in 
the present paper 



2.2. The background fluctuation distribution P(D) as a furtlier 
constraint on tiie faint source statistics 

The analysis of the fluctuations in the background intensity pro- 
vides other useful constraints on the faint source distribution in 
circumstances where the map's spatial resolution is poor, such 
that relatively few sources are directly detectable above the lim- 
its imposed by the instrumental and source confusion noises. 
With the term source confusion we refer to the photometric con- 
fusion in the terminology of Dole, Lagache, & Puget (2003). 
The distribution of background fluctuations in sky cells, often 
called distribution P(D) of the total flux (or "deflection" D) in 
the cell, includes information on the source number counts down 
to fluxes corresponding to one source per beam, well below the 
nominal confusion limit of at least 10-20 beams per source usu- 
ally adopted (Scheuer 1974). We defer the reader to a detailed 
account in Appendix A of the computational procedures to cal- 
culate the predicted P{D) based on a model for the differential 
number counts and under the assumption of a random distribu- 
tion of sources on the sky. 

An interesting case is offered by the Spitzer/MIPS obser- 
vations obtained by the FIDEL project (Frayer et al. 2009; 
Magnelli et al. 2009). The P{D) analysis offers an important 



consistency check of the counts in Fig. |2]at fluxes fainter than 
10 mJy, where the direct counts statistics is poor, and limited to 
the sources detected in GOODS-N by Frayer et al. (2006a). 

Figure |6] shows a comparison of our model distribution of 
fluctuations D with that of background deflections measured 
inside circular apertures of different radii (from 12 to 18 arc- 
sec) from the COSMOS 70 yum map (Frayer et al. 2009). The 
observed distribution of deflections D is considered here to 
stem from the convolution of the instrumental noise (assumed 
Gaussian) and the source confusion distribution. The lines cor- 
responding to the Gaussian noise distribution (rms of 0.62 
mJy/beam, blue curve) and the probability distribution of de- 
flections D due to source confusion are shown. The convolution 
of the two provides a good fit to the observational distribution. 
The instrumental noise of 0.62 mJy/beam corresponds to 0.033 
MJy/sr for our beams of 450 square arcsec, and the Gaussian in- 
strumental noise is essentially set by the "negative" part of the 
P(D) distribution (the side corresponding to the negative values 
of D) in Fig. |6] Our estimate of the source confusion noise of 
cr confusion ~ 0.63 mJy is larger than found by Frayer et al. (2009), 
0.35+0. 15 mJy, although only at the 2cr confidence level. We do 
not have a clear explanation for this difference, except perhaps in 
the Frayer et al. procedure to perform a simple quadratic subtrac- 
tion of the Gaussian instrumental noise from the total observed 
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Fig. 8. Galaxy LLP's at 15 /urn from Fang et al. (1998) and 
adapted from Xu et al. (1998, red open squares) in the low- 
luminosity regime. Lines correspond to our best-fit model de- 
scribed in Sect. 4. The contribution to the 15 fim LLP by type-I 
AGNs is shown as green long-dashed line, type-II AGNs plus 
moderate-luminosity starbursts make the green short-dashed line 
(type-II AGNs and starbursts are included in the same popula- 
tion on the assumption that in both classes the IR spectrum is 
dominated by starburst emission). The red line corresponds to 
the local LP of a population of high-luminosity starbursts domi- 
nating the source emissivity at high redshifts. The dotted line is 
the separate contribution of normal spirals, while the continuous 
line is the total LLP. 



P(D), instead of a full deconvolution as done here, more appro- 
priate for skewed distributions. 

This analysis of the fluctuation distribution is important to 
confirm the fast convergence of the counts fainter than a few 
mJy in the 70 /urn counts: a shallower convergence of the counts 
(see Pig.© would imply excess signal in the P(D) distribution at 
fi > 0.05 MJy/sr compared to what is observed in Pig.|6] The 
P(D) data reported in the figure impose quite a fast convergence 
of these counts fainter than 10 mJy. 

2.3. The CIRB background intensity 

Studies of galaxy evolution in the IR benefit by including data 
on the integrated emission of all cosmic sources, in the form of a 
diffuse background radiation (the CIRB) detected in COBP all- 
sky maps between 100 and 1000 fim (Puget et al. 1996; Hauser 
et al. 1998; Lagache et al. 2004). We defer to these papers and 
to Pranceschini et al. (2008) for reviews on the subject. 

A summary of the available data on the CIRB intensity /(v), 
compared with our best-fit model predictions, is given in Pig. 
121 Those in the far-IR and sub-mm are direct detections of the 
CIRB based on COBP maps, once the zodiacal and galactic 
foregrounds are subtracted. The mid-IR datapoints instead are 
based on simple integrations of the number counts in the respec- 
tive wavebands published by Plbaz et al. (2002), Papovich et al. 
(2004), and Pazio et al. (2004). 

These sub-mm CIRB detections are relevant for constrain- 
ing the high-redshift evolution of the galaxy IR emissivity. Two 
effects balance together here: that of the K-correction, which is 



particularly strong because of the very steep sub-mm spectra of 
all galaxy types (S [v] oc v^'^), increases the source contribution 
at the high-z. On the other hand, for a given source comoving 
emissivity pmiz) at redshift z, its contribution to the background 
intensity decreases proportionally to redshift as 

v/(v) = cp/R /47r(l +z). 



0.0<z<0. 200000003 
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Fig. 9. Galaxy LLP's at 24 /im from our analysis of the com- 
bined SWIRE-SDSS samples in the Lockman Hole (Vaccari et 
al. 2009). The contribution to the 24 yum LLP by type-I AGNs is 
shown as a green long-dashed line, type-II AGNs plus starbursts 
make the green short-dash line. Type-II AGNs and starbursts are 
included in the same population on the assumption that the IR 
spectrum is dominated in both classes by starburst emission. The 
dotted line is the separate contribution of normal spirals, while 
the continuous line is the total LLP. 



2.4. Tine iocai boundary: iuminosity functions 

The local luminosity functions (LLP's) for galaxies selected at 
different IR wavelengths set a fundamental constraint on the 
population's cosmological evolution. Various independent esti- 
mates of the galaxy luminosity functions at different IR wave- 
lengths have been reported. In the mid-IR, LLP's have been pub- 
lished by Rush et al. (1993), Xu et al. (1998) and Pang et al. 
(1998), all based on the 12 fim IRAS survey from which statis- 
tically complete and numerous all-sky samples have been ob- 
tained. However, as discussed in AP2001, these were affected 
by various uncertainties, including the IRAS photometric errors 
and the effect of local inhomogeneities, particularly those from 
the the local supercluster, implying an improper flux normaliza- 
tion and too steep a faint-end slope (e.g. inconsistent with the 
IRAS 60 iJin LLP by Saunders et al. 1990). These various ef- 
fects have been discussed by Pang et al. (1998), Xu et al. (1998), 
and AP2001: a corrected and re-calibrated 15 fim luminosity 
function based on these analyses is reported in Pig. [8] (see also 
AP2001). This re-estimate shows in particular a shallower low- 
luminosity slope that turned out to be consistent with the 60yum 
LLP by Saunders et al. (1990). 

We have obtained new improved estimates of the IR 
multi-wavelength luminosity functions by cross-correlating the 
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SWIRE survey (Lonsdale et al. 2003, 2004) in a large (10 
sq.deg.) area in the Lockman Hole with the Sloan Digital Sky 
Survey (SDSS, Stoughton et al. 2002) spectroscopic and pho- 
tometric redshifts. The combination of the high levels of pho- 
tometric accuracy, completeness, and reliability of the SWIRE 
source catalogues and the completeness of the SDSS redshift 
information provided us with a very reliable source list for lu- 
minosity function studies. Also, the large area covered by the 
combined surveys in Lockman appears free of conspicuous local 
galaxy condensations and offers an unbiassed sample of the local 
universe. While all details about this analysis will be reported in 
Vaccari et al. (in progress), we summarize here the main facts. 



10-' pT 



0.0<z<0.25 
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log L 

Fig. 10. Galaxy LLP's at 70 fj.m from our analysis of the com- 
bined SWIRE-SDSS samples in the Lockman Hole. The black 
filled square datapoints are as in Fig. |9] cyan open squares are 
the 60 /im local LP by Saunders et al. (1990) scaled to 70 fj.m 
with a typical galaxy SED. Lines as in Pig. [8] 

Vaccari et al. (20010, and in preparation) have considered 
all three SWIRE MIPS catalogues at 24, 70, and 160 fj.m and 
performed a completeness analysis of the three samples using 
Monte Carlo simulations. A likelihood ratio algorithm (Vaccari 
et al. 2008) was used to perform a cross-correlation of the 
SWIRE and SDSS samples. We also made use of the IRAC sur- 
vey in the whole field to eliminate spurious or uncertain identi- 
fications. The extension of the area covered by IRAC (9 sq.deg.) 
set the size of the combined survey. 

Eventually, the compromise MIPS flux limits have been set 
at 1.5, 20, and 90 mJy at 24/70/160 fim, respectively, values that 
ensure a high level of completeness. The redshift boundary for 
the computation of the local luminosity function was set in all 
cases to Zm = 0.2, that guarantees good statistics and dynamic 
range in luminosity and completeness in the redshift informa- 
tion. For z < Zm the fraction of spectroscopic redshifts is still 
large (> 50%) compared to the less secure photometric ones. 
With these boundaries, our final samples included 1058, 664, 
and 293 sources in the three MIPS channels at 24, 70, and 160 
fim, respectively. 

The local luminosity functions were calculated using the 
IfVmax estimator and Poisson error statistics. The K-corrections 
have been simply obtained with power-law fits to the average 
spectra in the narrow frequency interval defined by Zm, with flux 
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Fig. 11. Galaxy LLP's at 160 fim from our analysis of the com- 
bined SWIRE-SDSS samples in the Lockman Hole. Lines as in 
Fig. El 



densities assumed to scale as Sy oc v at 24 and 70 jim and 
constant with frequency at 160 yum. 

Preliminary results of this analysis are reported in Pigs. |9] 
[TOl and[TT] Our results in Pig. [TO] in particular, are compared 
with those obtained from the IRAS all-sky survey and later spec- 
troscopic follow-up by Sauders et al. (1990). This comparison 
shows a fairly good match of the two, except in our highest lumi- 
nosity bin, where our statistics are poor compared to IRAS. The 
same figures also report modellistic fits that will be discussed 
later 

2.5. The redshift distributions for faint IR sources 

Thanks to Spitzer's MIPS sensitivity, faint source catalogues 
from its surveys include substantial numbers of high and very 
high redshift objects. This is particularly the case for the 24 jum 
catalogues from GOODS, SWIRE, PES, COSMOS, and other 
surveys, while those at 70 and 160 fim achieve shallower depths 
and suffer from the more problematic optical identification. For 
the first time with Spitzer, we can infer relevant constraints on 
the galaxy and AGN infrared emissivity far above redshift of 
unity, including the main assembly phases of galaxies and of the 
cosmological structure. 

The simplest and most direct observational constraint on the 
high-redshift evolution can be inferred from the distribution of 
the redshifts for complete flux-limited source samples. The pro- 
cedure to obtain this simple information rests on the optical iden- 
tification of all sources in the sample and on a measurement of 
their (either spectroscopic or photometric) redshift. Once this is 
obtained, a comparison with the modellistic expectation for the 
redshift distribution is easily done via a simple integral of the 
redshift-dependent LP over luminosity. 

Rodighiero et al. (2010) have analysed multi-wavelength 
data, including deep photometric imaging with Spitzer, in three 
areas - GOODS-South, GOODS-North and SWIRE/VVDS - 
and have derived flux-limited samples at 24 yum. The sky areas 
are 168, 180, and 2160 square arcmin, and the limiting fluxes 
80, 80, and 400 juJy, respectively. The corresponding final num- 
ber of sources in the complete catalogues are 666, 889, and 
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Fig. 12. Left: The observed redshift distributions of 24 /vm sources with 5 24 > 80 jJLJy in the GOODS-S field (red histogram). Right: 
same as in the left panel but for a brighter flux limit 524 > 250 ^Jy in the total area of 348 sq.arcmin covered by GOODS-S and 
GOODS-N. 

Flim= 0.005000 Jy; Area (sqdeg)=0.200; lambda=850.0 mu 




S(850 micron) [mjy] 



Fig. 13. Statistics of SCUBA sub-mm sources. Left: The Euclidean-normalized diff'erential counts by Coppin et al. (2006) compared 
with the model expectation. Right: The observed redshift distributions of 850 |/m sources reported by Chapman et al. (2005, red 
histogram), compared with our model prediction for 5 850 > 5 mJy (arbitrary normalization). The datapoint at z=1.6 has been 
corrected to account for incompleteness in the spectroscopic survey. 



1495 for GOODS-S, GOODS-N and SWIRE/VVDS, of which 
64%, 33%, and 9% have spectroscopic redshifts. The remaining 
sources had reliable redshift estimates from photometric analy- 
ses, as explained in Rodighiero et al. These Spitzer 24 yum se- 
lected samples serve also as a basic dataset for the luminosity 
function studies in the next Section. 

The optical identifications of these 24 fxm source samples 
were performed on the basis of a maximum-likelihood technique 
taking the chance for spurious matches into account in detail. 
This procedure was particularly successful thanks to the avail- 
ability of very deep imaging with Spitzer/IRAC in the near-IR 



in all fields, providing a spectral bridge between the far-IR and 
the optical catalogues. As a result, only minor fractions of the 
Spitzer sources lacked an optical association. They were 3% and 
1.7% in GOODS-S and GOODS-N, while essentiafly afl sources 
in SWIRE/VVDS were reliably identified given the moderate 
depth of the primary IR selection. Fig. [T2] summarizes our re- 
sults for the redshift distributions in the deepest areas, compared 
with our model predictions, for two flux limits. 

On the side of sub-mm selected source samples, with their 
optical identification and spectroscopic characterization made so 
difficult by the optical faintness and large error boxes, the most 
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Fig. 14, Galaxy LP's in the rest-frame 15 //m in six redshift intervals. The blue thick hnes indicate the luminosity functions obtained 
by Rodighiero et al. (2010) in their combined GOODS+VVDS-SWIRE sample. The thick green dashed lines are the local LF by 
Pozzi et al. (2004) and are reported in all panels for comparison. The upper black dashed lines are our best-fit model prediction. 



systematic attempt to determine their redshift distribution is re- 
ported in Chapman et al. (2005) based on a combined SCUBA 
and very deep radio survey. The right-hand panel in Fig. [T3] 
shows a comparison of the Chapman et al. redshift distribution 
with our model's expectation. 



2.6. Direct estimates of the far-IR luminosity functions in 
redshift intervals 

Ideally, quite detailed statistical information on galaxy emissiv- 
ity and its cosmological evolution can be inferred from redshift 
surveys of complete flux-limited samples, like those described 
in the previous section. A popular procedure is to derive esti- 
mates of the galaxy luminosity functions in redshift intervals. 
In practice, this approach suffers various limitations. There is a 
first obvious need for very large source samples if one wants to 
achieve an acceptable resolution and to gather enough statistics 
in redshift space. More seriously, complete samples with (ei- 
ther spectroscopic or photometric) redshift information require 
not only an enormous effort of spectroscopic follow-up, but also 
ultra-deep imaging in complementary optical and near-IR wave- 
lengths to obtain reliable estimates of photometric redshifts for 
the - typically large - fraction of sources without any spectro- 
scopic detection. The implication is that a complete redshift cov- 
erage is only achieved over small sky areas and moderate depths, 
and appropriate sampling of the luminosity function is not possi- 
ble at low luminosity, except for very low-redshift bins. Finally, 



a representative sampling of the whole redshift-luminosity ob- 
servable plane would require combining ultradeep samples over 
small sky areas with shallower surveys on much larger fields (the 
so-called wedding-cake cosmological surveying). 

Attempts to reconstruct the galaxy and AGN far-IR lumi- 
nosity functions by various authors (e.g. Babbedge et al. 2007; 
Caputi et al. 2006, 2007) have exploited the extraordinary deep- 
mapping capabilities of Spitzer/MIPS at 24 yum. We made use 
of the results of a systematic analysis carried out by Rodighiero 
et al. (2010) on the basis of the optimal combination of depths 
and areas offered by the three Spitzer samples described in the 
previous section. 

While we defer to Rodighiero et al. for any detailed informa- 
tion, we report the results of that work in six redshift bins from 
z=0 to 2.5 in Figs.[T3and[T5] based on the combined GOODS-S, 
GOODS-N and SWIRE/VVDS sample including a grand-total 
of 3050 galaxies. Figs. [14] and [15] illustrate the data-to-model 
comparisons at 15 and 24 yum, respectively. In both figures the 
effects of the flux limits in the various surveys are apparent, pre- 
venting detections of faint sources at progressively higher red- 
shifts. If sources down to low luminosities (L ^ 10^ Lq) are 
detectable in the lowest redshift interval, this threshold luminos- 
ity increases to 10'" Lq by z=l and 10" Lq at z=2.5. We defer 
to Rodighiero et al. for an extensive comparison with other pub- 
lished results in the literature, showing a good overall agreement. 

Overall, there appears to be excellent agreement between the 
model and data, particularly ?A A^ff - 15 yum where the K- 
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Fig. 15. Galaxy LP's at 24 fim in six redshift intervals, (same as Fig. [141 except at 24 yum). In the low-redshft bin, a comparison is 
performed with results from Vaccari et al. (2009 in preparation), Marleau et al. (2007) and Shupe et al. (1998). The latter is used in 
all panels to show the effect of redshift evolution. 



correction for the data is minimal at the typical source redshift 
of z * 1 . The match is slightly less precise at A^/f - 24 yum, par- 
ticularly at high redshifts where K-corrections are substantial. 

Another aspect to note in the top-left panel of Fig. [TSj cor- 
responding to the local LF from the combined sample, is an ex- 
cellent match with the results of the SWIRE/SDSS analysis in 
Sect. 12.41 and Fig. |9l The two datasets seem to be in slight ex- 
cess compared with the model, perhaps due to some large-scale 
inhomogeneities in the small local volume sampled. 

We also compare in Fig.[T6]our model's results with the lu- 
minosity functions in redshift bins based on the FIDEL sample 
by Magnelli et al. (2009). This sample is based on very deep 
maps with Spitzer/MIPS in GOODS-N, GOODS-S and ECDS 
and includes 680 sources selected at 70 fim with almost com- 
plete photometric or spectroscopic redshift information. As ap- 
parent in the figure, excellent agreement with the model is found 
in this case too. 



3. Multi-wavelength modelling of IR galaxy 
evolution 

Our adopted multi-wavelength modelling of the IR galaxy 
evolution follows a similar approach to the one pursued by 
Franceschini et al. (2001, AF2001), which was calibrated on a 
much shallower database essentially restricted to the ISO sur- 
vey data. Compared with that scheme, some significant improve- 
ments were required, as discussed later. 



The present multi-wavelength analysis is not intended to pro- 
vide us with a detailed and complete justification of the complex 
physical processes involved in star-formation and galaxy assem- 
bly and traced by the IR emission. Instead the model guarantees 
the best possible adherence to an enormous variety of data on the 
cosmic source emission from the local to the distant universe and 
is intended to identify general patterns in galaxy evolution use- 
ful for comparison with more physically -based galaxy formation 
models. 

3.1. The population composition 

Our model population includes contributions by various source 
categories characterized by different physical and evolutionary 
properties. 

1 . The first considered class of sources are non-evolving nor- 
mal spirals that dominate the multi- wavelength LLF's at low 
luminosities. Such a population was required in the AF2001 
analysis to explain some features in the extragalactic source 
number counts, particularly those at 15;um available at that 
time: an Euclidean slope of the counts down to the milli- 
Jy level, with a sudden super-Euclidean excess appearing at 
fainter fluxes. This behaviour is completely confirmed by the 
new Spitzer/MIPS 24/im counts shown in Fig. [1] and is also 
apparent in the Spitzer data at 160 /vm and partly those at 
70 jum in Figs.[3]and|2] These data are inconsistent with the 
whole local galaxy population to evolve back in cosmic time 
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Fig. 16. Galaxy LP's at 35 jum in three redshift bins, based on the FIDEL 70 jum selected sample by Magnelli et al. (2009). 



and instead they require the summed contributions of two 
source categories, a non-evolving and a fast-evolving popu- 
lation. We identify this non-evolving component as a popula- 
tion of spiral galaxies, building up progressively with cosmic 
time, as indicated by the rather flat age distributions of their 
stellar populations observed in local objects. We assumed for 
this class of sources a redshift of formation of z form - 1 and a 
local normalization consistent with the observed local lumi- 
nosity functions discussed in Sect. 2. Our adopted IR spectral 
energy distributions (SEDs) for spiral galaxies are discussed 
in Sect. [121 

2. A second obvious class of sources to be considered in the 
model and characterized by high rates of evolution are active 
galactic nuclei, which are found in substantial numbers in 
IR selected surveys at all wavelengths (Franceschini et al. 
2005; Rowan-Robinson et al. 2005). Among these, type-I 
AGNs are straightforward to identify based on simple com- 
binations of optical-IR colours (Hatziminaoglou et al. 2005; 
Siana et al. 2008). For the complementary category of type- 
II objects instead, the identification would at least require the 
combination of far-IR/sub-mm data with very deep optical, 
near- and mid-IR imaging to detect the very hot circumnu- 
clear dust emission characterizing them. For this and other 
reasons, this class of sources will be treated separately later 
The local luminosity function for type-I AGNs has been 
taken from AF2001, which was essentially fitting data by 
Rush et al. (1995) based on the IRAS 12 yum survey. Our as- 
sumed LLF, shown as green dot-dashed line in Fig.|8]and the 
following, is entirely consistent with the Rush et al. Seyfert- 



1 LLF. The LLFs at other IR wavelengths are also reported 
inFigs.milOlIII] 

Type-I AGNs and quasars are modelled based on sim- 
ple evolutionary recipes inferred from optical and X-ray 
quasar surveys and with standard UV-optical-IR-mm spec- 
tra. As a guideline for assessing their evolution properties in 
the infrared, we have used the catalogue of 1493 SWIRE 
sources selected at 24 yum in an area of 0.85 sq. degree 
in SWIRE/VVDS to a flux limit of 5 24 = 400 A^Jy as re- 
ported by Rodighiero et al. (2010, see Sect. 12.6b . A spectral 
analysis making use of the optical and IR data available in 
the field for all objects in the complete sample is reported 
by Rodighiero et al., among which the type-I AGN (that 
is quasars and Seyfert-1 galaxies) are easily identified by 
their flat spectral shapes over the optical through IR wave- 
length interval. About 8.3% of the sources at the flux limit 
are found to be type-I QSO or Seyferts. We have used the 
24 yum number counts from this catalogue to constrain the 
evolution properties of Type-I AGNs in our model. In addi- 
tion, the bright ends (where AGN contribute significantly) 
of the redshift-dependent luminosity functions at 24 and 35 
yum discussed in Sect. 12.61 are used to further constrain this 
evolution. Based on these data, type-I AGNs are assumed to 
evolve in luminosity as L(z) - L(Q) x (1 -H z)^'^ up to z = 2, 
with L(z) - L{0)x3^'^ above. At the same time, these data in- 
dicate for the population a moderate comoving density evo- 
lution like dN/dL oc (1 + z)°-^ up to z=2. 
3. As discussed by AF2001 and Fadda et al. (2002) among oth- 
ers, AGNs of any kinds are far insuflicient to explain the 



Franceschini, A., et al.: Galaxy Evolution at IR Wavelengths 



13 



faint IR source statistics, because their high IR luminosities 
contribute to the counts only at rather bright fluxes (S 24 > 
10 mJy), while they are too rare to explain the large observed 
numbers at fainter limits. To try explaining the fainter num- 
ber counts we need to refer to a much more numerous pop- 
ulation of fast evolving objects. A third component was then 
considered, starburst galaxies of moderate luminosities, or 
LIRGs, which are already present in AF2001, Chary & Elbaz 
(2001), Elbaz et al. (2002), Bales et al. (1999), Lagache et al. 
(2003) and Takeuchi et (2006), among others. The increas- 
ing relevance of this population with redshift has been one 
of the main outcomes of the ISO mission and has been con- 
firmed by the Spitzer/MIPS observations (Le Floc'h et al. 
2005; Perez-Gonzales et al. 2005). All these analyses found 
that the MIPS data require a fast increase in the comoving 
emissivity up to at least z - I- 

The local fraction of the evolving moderately-luminous star- 
burst population is assumed to be ~ 10 percent of the total 
galaxy population (see Fig. [8]l, roughly consistent with the 
local observed fraction of interacting galaxies. The quick ob- 
served upturns in the normalized differential number counts 
in Figs. [T]to|4]then require a strong increase in the average 
emissivity of the evolving population, an increase with red- 
shift that is obtained with both luminosity and number den- 
sity evolution: 



dN , 



dL 



dN 



-Tr(Ll2[z],z) = \ — \ (Li2)x(l+Zi) 



dL 



•,5.2 



Z <Zl 



Zl < z < z„ 



Li2(z) = Ll2X(l+z)2'' z<Z2 



Lniz) = L12 x(l+ zif-'^ Z2<z<Zi 



(1) 



where ^(Li2[z],z) is the redshift-dependent luminosity 
function at our reference wavelength of 12 //m, and where 
Z\ - 0.72, Z2 - 0.63 and z-i - 1.2. For z > Zi the average 
population luminosity converges like L^iz) °^ (1 +zY^, up to 
Zmax - 5. This detailed pattern for the evolution of moderate- 
luminosity starbursts will be discussed later in comparison 
with that of higher-luminosity sources. 
4. Current data impose an important modification to the 
AF2001 evolutionary scheme: they require the introduction 
of a fourth evolutionary population component of very lu- 
minous starburst galaxies, or ultra-luminous infrared galax- 
ies (ULIRGs), dominating the cosmic IR emissivity above 
z - 1.5. This latter emerges in particular from the analysis 
of the galaxy samples selected by the deep Spitzer/MIPS 24 
surveys, and their improved capability of detecting cos- 
mic far-IR sources at redshifts z > 1 and up to z ^ 3 (see 
Fig. [TTlbelow). As discussed later in Sect.|4] these observa- 
tions have revealed the existence of a separate population of 
ULIRGs at high-redshift, essentially absent or very rare lo- 
cally, hence characterized by an extremely fast evolution in 
cosmic time. Millimetric observations with SCUBA/JCMT 
since several years, confirmed by other more recent observa- 
tions, have also found clear evidence of such very luminous 
objects at high redshifts (see Figs. [T2] and [T3ll. 
In our interpretative scheme the distinction between LIRGs 
and ULIRGs is not to be considered as literally related with 
their luminosities at any redshifts. Indeed, as shown in Fig. 
[19] below, at low redshifts their luminosity distributions are 



similar, but they become quite distinct at z > 1, where LIRGs 
have peak luminosity around Lt,,,! - 10^' Lq, and ULIRGs 
around Lboi - 10'^ L©. 

Our modelling of this population assumes that it makes a 
negligible contribution to the local LF's (shown as red long 
dash lines in Figs.lSltofTTTi. mainly confined to the highl lumi- 
nosity end of the LLF's, and then evolve in time in a similar 
fashion as eq.[T] 



Z < Zl 



dN ldN\ 

^(ii2[z],z) = ^^j (Li2)x(l+zi)4-> zi<z<z„ 



Ln{z)^Ln>^{l+zt^ z < Z2 



il2(z) = il2 X (1 H-Z2)^'' Z2<Z<Z3 



(2) 



where ^(ii2[z],z) is the redshift-dependent ULIRG lumi- 
nosity function at our reference wavelength of 12 yum, and 
where zi = 1 -3, Z2 — 1-3 and Z3 - 1.2. For z > Z3 the average 
population luminosity converges like Li2(z) oc (1 H- z)"^^. 

Our assumption here is to include in the two starburst cat- 
egories also type-II AGNs of moderate or high IR luminosity 
(type-II Seyfert galaxies and dusty quasars), whose properties 
of nuclear emission are always associated with various levels of 
circum-nuclear starburst activity. Disentangling the two physical 
emissions, particularly in faint distant sources, is a major prob- 
lem of current-day cosmology (see e.g. Daddi et al. 2007 for 
Spitzer sources and Alexander et al. 2005 for SCUBA objects). 
In addition to a close relationship of starburst and dust-obscured 
AGN emission observed in type-II AGNs, there is also evidence 
of a similar rate of cosmological evolution at least up to redshift 
z=l for the two source categories (e.g. Franceschini et al. 1999). 
Finally, the IR spectra of starbursts and obscured AGNs reveal 
rather similar shapes (e.g. Berta et al. 2003; Polletta et al. 2007), 
while diff'erences amount to a slight excess mid-IR emission in 
dusty AGNs. For all these reasons, and in the spirit of keeping 
the model complexity and number of free parameters to a mini- 
mum, we included type-II AGNs into the starburst categories. 

The AGN contribution to the LIRG and ULIRG classes is a 
major question discussed in Sects. l4~2l and l4.4l below. 

3.2. The spectral model 



As anticipated in Sect. 12.41 and following AF2001, our multi- 
wavelength modelling scheme assumes as a basic reference the 
local luminosity functions of various populations based on the 
IRAS 12 jum all-sky survey. This choice is justified by the fact 
that the IRAS survey was, and still is, the only all-sky IR sur- 
vey for which a complete source identification and spectroscopic 
follow-up has been obtained, allowing us a reliable estimate of 
the galaxy local IR luminosity functions. Fig. [8] shows a com- 
parison of our modellistic LLF and the observed datapoints at 
the closeby eff'ective wavelength of 15 fim, showing good agree- 
ment of the two. From this LLF, luminosity functions at all IR 
wavelengths are then computed based on a broad-band spectral 
model. As for the latter, our choice was again to keep to a mini- 
mum the complexity of the model. 

In consideration of the multi-wavelength nature of our used 
database, the definition of a spectral model is particularly im- 
portant. Our choice has been to define a-priory only the average 
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Fig. 17. Our adopted IR spectra of various galaxy populations. The short-dashed cyan line corresponds to our adopted spectrum for 
the moderate-luminosity LIRG starburst population, while the red long-dashed curve is the spectrum of high-luminosity ULIRG 
sources. In both cases the spectra are similar to that of the prototype star-forming galaxy M82 (in the range from 5 to 18 jum it is 
precisely the ISOCAM CVF spectrum of M82). The lower dotted line corresponds to a low-luminosity inactive spiral {Siov.\y\ in 
eq. IA.2I ). while the upper dotted line is closer to that of ULIRGs {S higi, [v]). The lower dot-dash green line is the average type-I AGN 
spectrum. The boundaries of the MIPS 24 yum filter are also shown in the source rest-frames at various source redshifts: due to the 
prominent PAH features, we expect relative maxima in the redshift distributions at z ~ 1 and ~ 2, in agreement with the observations 
(Fig. [13. 



spectrum of the type-I AGNs, because this has little influence 
on the model outcomes except for the 24 fim statistics, as pre- 
viously discussed in Sect. 13.11 For the other population compo- 
nents (spirals, LIRGs and ULIRGs, with the exception of the 
low-luminosity spirals which have cooler spectra), we have ba- 
sically referred our model spectra to that of the prototype local 
starburst galaxy, M82. We have then made adjustments to it so 
as to achieve a best-fitting to the available data from the mid-IR 
throughout the millimetre. We anticipate that an interesting out- 
come of our analysis was that small modifications to the M82 
broad-band spectrum are sufficient for obtaining this best-fit. 

3.2.1. AGN 

For the type-I AGN population we assume an SED correspond- 
ing to an emission model by a face-on dusty torus and is reported 
as a green dot-dashed line in Fig. [17] This average spectrum is 
taken from the radiative transfer model by Fritz et al. (2006), and 
was optimized to fit photometric data for AGNs in the SWIRE 
surveys (see Franceschini et al. 2005; Hatziminaoglou et al. 
2008). 

Our adopted type-I AGN spectrum is close to and consistent 
with the average Spitzer/IRS spectral shapes of quasars pub- 
lished by Netzer et al. (2007). A discussion is reported there 
about the exact origin of the far-IR flux of quasars and an in- 
dication is provided that this long-wavelength part, as opposed 
to the mid-IR spectrum which is certainly due to AGN circum- 



nuclear dust, might instead be due to stellar-illuminated dust. For 
our present purposes it makes no much difference if part of the 
AGN population radiation is due to starburst rather than gravita- 
tional processes. Note finally that our detailed spectral modelling 
for AGNs is not critical, because type-I AGNs contribute signif- 
icantly only to the 24 /im, and to a lower extent to the 70 jum, 
statistics at the bright fluxes. 



3.2.2. Spirals 

Also in the case of the normal spiral population a precise defini- 
tion of the spectral model is not critical for our attempt to inter- 
prete the high-redshift source evolutionary properties. Although 
spirals dominate the low-luminosity end of the LLFs, their lack 
of evolution implies that their number density becomes negli- 
gible compared to other evolving populations at cosmological 
redshifts. As discussed in AF2001, a comparison of LLF's in 
the mid-IR, far-IR and sub-mm shows that the 60yum LLF has a 
flatter power-law shape at high-luminosities compared with the 
other LLF's, an effect of the spectra for more luminous galaxies 
showing excess 60 //m emission compared with inactive galaxies 
(see also the luminosity-dependence of the IRAS far-IR colours 
of galaxies; Soifer et al. 1989; Chary & Elbaz 2001). We have 
then modelled the redshift-dependent multi-wavelength LLF's 
of spiral galaxies by assuming SEDs dependent on luminosity, 
with spectra ranging from those typical of low-luminosity inac- 
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live objects with peak emission at 200-300 jum, to those peaking 
at 60-80 fim characteristic of luminous luminous IR galaxies. 

We took as reference for our multi-wavelength LF the one 
by IRAS at 12 jjin (^{Ln)) discussed in Sect. 5.1. This is 
transformed to longer wavelengths according to SEDs that vary 
according to the value of L12. The assumption was that for 
L12 < 2 10^ Lq the spectrum is that of an inactive spiral (5/o,v[v], 
the lower black dotted line in Fig.fTTI). while for L12 > 10'^ Lq 
it is a typical ULIRG spectrum {S higii[v], upper dotted line in 
Fig.fTTli. For intermediate luminosity objects, the assumed SED, 
S (v), is a linear interpolations between the two as a function of 
the L12 luminosity: 

„ 5/o,,(v)(logLi2-9.3) + 5/,/,;/,(v)(logLi2-12) 
S (v)Il,2 = • (3) 

The multi-wavelength LLFs at any wavelengths A are then easily 
computed from the relation 

—{La\Li2,z) ^ —(Lx2,z)\— —\ . (4) 

dL dL \d10gL12l 

The 5(v)|lp spectrum from eq.[3]is also needed to compute the 
K-correction. 

3.2.3. LIRG and ULIRG galaxies 

Of course, more critical are the assumptions about the spectral 
fitting of the two luminous galaxy populations considered in the 
model, the moderate and the high-luminosity objects (i.e. the 
LIRGs and ULIRGs). To simplify our treatment, we adopted for 
each one of the two classes a single spectral form, independent of 
luminosity. For the LIRG and ULIRG objects our best-fit spectra 
are shown in Fig. [17] It should be noticed how similar the two 
spectra are, both of them turning out to be close to the spectrum 
of the prototype local starburst galaxy M82 (e.g. Polletta et al. 
2007). 

We then support earlier findings by AF2001 that an M82- 
like IR SED is needed to reconcile the various statistics of faint 
starbursting sources, including those currently detected in large 
numbers at z > 1 not sampled by previous ISO observations. 
In particular, the adoption of the IR SED of the ULIRG galaxy 
Arp 220, an highly absorbed ultra-luminous source, would im- 
ply the far-IR counts and EE's having far too high normaliza- 
tions compared with the mid-IR ones (at 12, 15, and 24 fim), 
and a complete mismatch of the observational data. The aver- 
age spectral shapes in Fig. [17] are required to reconcile a vast 
amount of multi-wavelength data and should then be considered 
as representative of the bulk of the population. 

We stress that, in spite of the simplifying assumption about 
the average spectra of IR-selected starbursts to be independent 
of redshift and luminosity, our model is still able to provide 
good fits to the various statistical observables. Obviously this 
is not by itself proof that the IR spectra are completely time- and 
luminosity-invariant. 

3.3. Model to data confrontation 

With our spectral model for the different source populations, we 
computed LFs at all IR wavelengths by convolving the spectra 
with the detailed filter response functions of the various photo- 
metric channels. Once the redshift-dependent multi-wavelength 
luminosity functions for our source populations are set, the pro- 
cedures for comparing them to the variety of data discussed 
in Sect. 2 are also straightforward. We defer the reader to a 



detailed account in AF2001 about the computation of number 
counts, redshift distributions, etc. Here again our spectral model 
is needed in order to calculate not only the luminosity func- 
tions at various wavelengths, but also the K-corrections in the 
luminosity-distance relation. All this requires that the average 
source spectra are convolved with the detailed spectral response 
functions of the various observational filters, as discussed in 
AF2001. 

As illustrated in Figs. \T\ to [16] our best-fit model appears 
to reproduce a wide variety of multi-wavelength IR data, from 
redshift-dependent luminosity functions to faint number counts, 
and providing marginally acceptable fits of the background spec- 
tral intensity and cell-to-cell fluctuations. A detailed comparison 
of our model for IR galaxy evolution with the luminosity func- 
tion data in redshift bins at 15, 24, and 35 fim (Sect. l276l l shows 
excellent agreement up to the maximum observed redshift of 
z - 2.5, with overall acceptable values of the total reduced 
(xl - 1 ■ !)■ There is only a slight tendency of the model fit to re- 
main low in the lowest bin at 24 yum and < z < 0.25, probably 
due to the unaccounted effects of the local large-scale structure 
in the luminosity function data. In any case, the low-z behaviour 
of the model was calibrated on a variety of number count data 
at bright fluxes (including the IRAS all-sky far-IR counts and 
those from the Spitzer SWIRE project), which guarantees excel- 
lent control of the local IR universal emissivity. 

Compared with published analyses of IR galaxy evolution 
based on direct estimates of the galaxy luminosity functions (e.g. 
Perez-Gonzales et al. 2005; Caputi et al. 2007; Rodighiero et al. 
2010), our current approach of combining these with integrated 
information coming from faint number count and background 
data allows us to better constrain the luminosity functions over 
wider redshift and luminosity intervals, including some extrapo- 
lations to the low luminosities where they are not directly mea- 
surable (see Fig. I20lbelow for more details about uncertainties 
in this analysis). 

3.4. Model predictions 

Appendix B is dedicated to discussing predictions for number 
counts by forthcoming observations, particularly with Herschel. 
We also report there some expectations for source confusion lim- 
itations with future instrumentation. 

Fits to the observational data and individual population com- 
positions are reported in Figs.[T]through[T8]with reference to our 
best-fit model. 



4. Discussion 

4.1. Generai evoiutionary patterns for IR sources: evidence 
of a strongly luminosity-dependent evolution 

The redshift evolution of the comoving bolometric IR emissivi- 
ties from 6 to 1000 /im for our model is reported in Figure [T8] 
Individual contributions of our four galaxy classes, the normal 
spirals, type-I AGNs, moderate-luminosity, and high-luminosity 
galaxies are also shown. Our model predictions are compared 
with a recent estimate of the IR emission of the global population 
based on 24 //m data and on fitting the LF data in redshift bins 
with two-power-law functions (Rodighiero et al. 2010). Both our 
model and the data of Rodighiero et al. show a fast increase in 
the comoving IR emissivity to z ~ 1, the total increase for the 
best-fitting model amounting to a factor ~ 30. A flattening of 
the comoving emissivity is then observed thereafter, while some 
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Fig. 19. The comoving bolometric emissivity from 6 to 1000 yum as a function of luminosity for the IR-selected galaxy population, 
calculated at various cosmic epochs from z=0 to 2.5. The IR emissivity is expressed here in solar luminosities per cubic Mpc per 
unit logarithmic interval of L. Line types as in Fig.fTSl 



convergence would be expected by the model at z > 2.5, as dis- 
cussed below. 

Figure [18] also reveals some remarkable differences in the 
evolutionary properties of the two energetically-dominant popu- 
lations of moderate- (LIRG) and high-luminosity (ULIRG) ac- 
tive galaxies: the epoch of peak activity of the former appears 
to be shifted to lower redshifts compared to the more luminous 
ones dominating at the highest redshifts. Similar findings based 
on IR source populations are also reported by Le Floch et al. 
(2005) and Perez-Gonzales et al. (2005). 

This general evolutionary pattern is confirmed by the more 
detailed analysis reported in Fig. [191 showing the evolution of 
the comoving bolometric emissivities from 6 to 1000 yum as a 
function of luminosity for the various IR-selected galaxy popu- 
lations, calculated in six redshift bins from z=0 to 2.5. The IR 
emissivities on the y-axis are expressed here in solar luminosi- 
ties per cubic Mpc and per unit logarithmic interval of L^oi, and 
are simply given by the product of the bolometric luminosity 
function times Lj,,,! (hence proportional to the contribution of the 
population to the total emissivity per unit decade of L). 

In the lowest-z bin, the bolometric emission is dominated by 
inactive spirals, type-1 AGNs contribute at the highest-L, while 
the ultra-luminous starburst ULIRG population provides a negli- 
gible contribution at any Lhoi- This situation rapidly changes al- 
ready by z ^ 0.5, when active galaxies of both LIRG and ULIRG 
classes start to dominate, except at the lowest values of L. The 
moderately-luminous objects (LlRGs) then appear to stop evolv- 



ing by z ~ 1, whereas the ULlRGs continue to rise in both co- 
moving number and bolometric luminosity to higher z. 

Remarkably enough, the highest redshift bins show a promi- 
nent dominance of the ULIRGs, with rather steady comoving 
number densities as a function of redshift. The LIRG population 
keeps a bumpy luminosity function at all redshifts, and already 
appears to decrease in number density and luminosity at z > 1 . 

A different view emphasizing the redshift-dependence of the 
evolution is reported in Fig. [20] showing the comoving number 
density of galaxies versus redshift for various luminosity classes, 
indicated on the right hand side of the figure (these are the aver- 
age Lboi values in the bin). The figure also includes a detailed as- 
sessment of how the main datasets constrain our solution in vari- 
ous luminosity and redshift intervals, thereby providing informa- 
tion about the reliability of our solutions. Continuous lines are 
regimes where the luminosity functions are constrained by cur- 
rent far-lR number counts, particularly the Spitzer 70 jum counts 
at the current sensitivity limit, in addition to all other avail- 
able statistics; the dashed lines are constrained by the Spitzer 
24 yum counts at the current sensitivity limit (30 yuJy), 24 //m 
luminosity functions and z-distributions. Inside these redshift 
intervals spectral extrapolation to the far-lR is required to esti- 
mate the bolometric luminosity, but this is also well constrained 
by BLAST and preliminary Herschel observations (Sect. lTTTSl l. 
Finally, dotted lines are regimes where the IR luminosity func- 
tions are constrained essentially by the CIRB spectral intensity 
and background fluctuations. As the available constraints de- 
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Fig. 20. The best-fit model prediction about the comoving number density versus redshift of cosmic IR sources as a function of 
source luminosity. The latter, indicated at the right hand side of each curves, is given by vL(v) at 100 fim in solar luminosities. Line 
types correspond to redshift intervals with the following coding. Continuous black lines are regimes where the luminosity functions 
are constrained by current far-IR number counts, particularly the Spitzer 70 yum counts at the current sensitivity limit, in addition to 
all other available statistics; the red dashed lines are constrained by the Spitzer 24 yum counts at the current sensitivity limit (30 /iJy), 
24 yum luminosity functions and z-distributions; cyan dotted lines are regimes where the IR luminosity functions are constrained 
essentially by the CIRB spectral intensity. As the available constraints decrease going from low- to high-z, the uncertainties in the 
curves for the various luminosity intervals also increase with redshift. 



crease going from low- to high-z, the uncertainties in the curves 
for the various luminosity intervals also increase with redshift, 
mostly because of an increase in the degeneracy of the model 
solutions. 

Again interesting trends are shown in Fig. |20l The rate of 
increase from z = to 1 of the comoving density is a strong 
and monotonic function of luminosity, the rise factor being more 
than 3 orders of magnitude for the most luminous galaxies, while 
it is only a factor ~ 3 for the faintest ones. (The previously men- 
tioned average factor for the total population is ~ 30). The other 
main observed feature is a global and systematic shift of the peak 
emissivity as a function of redshift, the most luminous objects 
having their maximal activity at z > 2, while the least luminous 
ones do it at redshift close to 1 . This details more precisely in 
terms of the bolometric luminosity class a trend already antici- 
pated in Fig.fTSl 

The lowest luminosity class in Fig.l20l(contributed bv spirals 
and LIRGs) reveals a moderate evolution over the whole redshift 
range analysed here. It should be noted, however, that the high- 
redshift behaviour of sources with this luminosity is not reliably 
sampled by current observations. For all other luminosity bins 
with Lboi > 10"^ Lq, instead, the combined dataset of luminosity 



functions, deep number counts, and background data considered 
by our analysis provides a relevant set of constraints on the evo- 
lution of emissivity of IR galaxies that are not available from 
galaxy surveys at any other wavelengths. 

Finally, the convergence of the comoving emissivity at z > 
2.5 in Fig. [18] was imposed in particular to avoid exceeding the 
SCUBA number counts, particularly the two faintest datapoints 
at 5 850 between 2 and 6 mJy in Fig.[T3]left, and to a less extent 
the CIRB COBE-observed intensity at /I > 500 yum. The function 
piR is dominated here by the ULIRG population. 

The detailed behaviour of the high-z convergence of the 
ULIRG's comoving volume emissivity (Fig.fTSl) is not yet well 
constrained by the scanty information on the z-distribution of 
SCUBA sources (Fig. [HI Chapman et al. 2005). Equally, the 
deep Spitzer/MIPS 24 yum surveys suffer a strong K-correction 
penalty at z > 3 (Fig.[T2li. such that a few Spitzer sources are de- 
tected out there. Confirmation, particularly by SCUBA-2, of the 
source number counts at the faintest levels, and the characteriza- 
tion and redshift distributions of the forthcoming large samples 
of Herschel sources in the sub-mm, will be essential to establish 
the detailed cosmic evolution at z > 2.5. Also later high spatial- 
resolution studies with ALMA will offer us richer datasets to 
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Fig. 18. Evolution of the comoving bolometric luminosity den- 
sity from 6 to 1000 jum for the IR-selected galaxy population, 
based on the model discussed in Sect. |3] The luminosity den- 
sity is expressed here in solar luminosities per cubic Mpc. Green 
dot-dashed line: type-1 AGNs. Cyan short-dashed line: evolving 
moderate-luminosity starbursts (LIRG). Red long-dashed line: 
high-luminosity starbursts (ULIRG). Lower dotted black line: 
quiescent spiral population (note that the spirals are assumed not 
to evolve in comoving luminosity by our model and then to dis- 
appear at z-1). The upper continuous line is the total predicted 
emissivity. Red filled square datapoints are from the analysis of 
the 24 fim luminosity functions by Rodighiero et al. (2010). Blue 
datapoints are from a large spectroscopic survey at 1.9 < z < 3.4 
by Reddy et al. (2008). 



constrain this early history of galaxy IR emissivity. We will as- 
sume henceforth that piR quickly converges there, as in Fig.fTSl 
but will later discuss the uncertainties related with this. 

Even when taking this uncertainty into account, there is a 
clear tendency in our redshift-dependent bolometric comoving 
energy density shown in Fig.[T8]to stay lower than various pub- 
lished results, particularly at higher redshifts than z ~ 1.5. These 
include Hopkins (2004), Hopkins and Beacom (2006), and the 
overview reported in Perez-Gonzalez et al. (2008), and are typ- 
ically based on either large extinction corrections (for the UV- 
optical, or Ha observations) or large spectral extrapolations (like 
those making use of the MlPS/24//m, the SCUBA, or deep radio 
surveys). This difference will have important implications for 
our later investigation. 

In conclusion, the evolution of IR emissivity for our best- 
solution model emphasizes the activity of very luminous objects 
at the higher redshifts (z ~ 2.5), while lower luminosity galax- 
ies appear to dominate at lower redshifts. This kind of evolution 
is very reminiscent of the downsizing scenario that has emerged 
from studies of the redshift-dependent galaxy stellar mass func- 
tions (Bundy et al. 2006; Franceschini et al. 2006; Fontana et 
al. 2006; Perez-Gonzalez et al. 2009; among others), in which 
the most massive galaxies and their stellar content are already 
largely in place at z > 1, while less massive systems continue to 
form stars until more recent epochs. We seem to have identified 
a corresponding evolutionary feature in terms of the e.m. (bolo- 



metric) emission produced by forming stars (and AGNs). This 
also confirms and quantifies evolutionary trends anticipated by 
Cowie et al. (1996) for galaxies and Franceschini et al. (1999) 
for both galaxies and AGN (see also Perez-Gonzalez et al. 2008; 
Mobasher et al. 2009). 

However, in terms of bolometric, comoving energy density, 
this activity appears to decrease at z > 1.5 and to remain slightly 
lower than sometimes reported earlier, based on large extrapola- 
tions. 

4.2. Origin of the IR emission in the sample sources 

The big advantage of far-lR and sub-mm observations of faint 
high-redshift objects is that, particularly for the more luminous 
source populations, the LlRGs and ULlRGs, these data offer an 
excellent proxy for the bolometric emission (e.g. Soifer et al. 
1987; Sanders and Mirabel 1996; Elbaz et al. 2002; Rieke et al. 
2009). If we know whether the emission is primarily to ascribe to 
either star-formation or AGN activity, we can then directly infer 
the evolutionary history of star formation or nuclear gravitational 
accretion from the IR emissivities. 

The problem here is that various different processes tend to 
mix together in producing IR photons, so that the interpretation 
may not be trivial. Fortunately, good insight into the physical 
nature of the primary energy sources in both local and high- 
redshift IR objects have been obtained with the very sensitive 
IRS infrared spectrograph on Spitzer. Yan et al. (2007), Sajina et 
al. (2007, 2008), and Lacy et al. (2007), among others, have re- 
ported low-resolution spectroscopy in the critical region between 
5 and 40 fim, where the stellar and AGN components show 
maximal differences. The IRS spectroscopy of high-redshift IR 
sources has been possible for the highest luminosity objects, 
which tend to include a higher incidence of AGNs (Vieilleux 
etal. 1997; Kim etal. 1998). 

As a general trend, these studies indicate that high-z sources 
selected in the sub-mm (the SCUBA galaxies) are dominated by 
starburst emission, while those selected by Spitzer/MIPS in the 
mid-IR tend to include a higher incidence of AGNs. The rela- 
tively higher contribution to the mid-IR than to the sub-mm by 
very hot circum nuclear dust heated by a nuclear AGN is well 
established by observational evidence and theoretical expecta- 
tion. Although including AGN contributions in a large major- 
ity of cases, the extensively investigated SCUBA galaxies ap- 
pear to have their bolometric flux dominated by star formation 
(Alexander et al. 2005). Daddi et al. (2007) from the analysis 
of luminous 24 pm selected galaxies at z ~ 2 find that the ma- 
jority of the sources more massive than 10'^ Mq show excess 
hard X-ray emission in stacked Chandra images, revealing trace 
emission by a Compton-thick obscured AGN. In general, how- 
ever, all X-ray studies of the AGN incidence in obscured high- 
redshift objects have to be taken with caution, due to the funda- 
mental limitation that, above some value of the column density 
(Nh - 10^^ cm^^), the whole X-ray flux is essentially repro- 
cessed to longer wavelengths. In this regard. X-ray studies pro- 
vide lower limits to the AGN content. 

An interesting analysis of the AGN incidence in IR- and sub- 
mm-selected sources has recently been reported by Nardini et al. 
(2008) for local ULIRGs and Watabe et al. (2009) for high-z ob- 
jects. Their conclusions, based on stacked Spitzer/IRS spectra 
of samples of sources selected in redshift and luminosity bins, 
may be summarized as follows. The majority of both sub-mm 
and 24 pm sources show evidence of some AGN contribution 
in the IRS 6 to 8 pm rest-frame spectra, this fraction becoming 
dominant for the 24 pm selected objects. In terms of the bolo- 
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metric IR-mm luminosity, the AGN contribution appears to mix 
with what is coming from a starburst. In SCUBA galaxies the 
overall AGN contribution is estimated at about 7%, while for 
the ultra-luminous 24 fim sources it settles around one third of 
the bolometric source emission, fairly independently of redshift. 
It should be kept in mind that all these results are to be taken as 
very preliminary: not only are they based on a handful of objects, 
but the spectral analysis was also performed on the stacked IRS 
spectra. 

Our evolutionary scheme includes a limited ability to dis- 
entangle the main different primary power sources, young stars 
in and around molecular clouds, and AGN emission from grav- 
itational accretion. The type-I broad-lined AGNs are in general 
easily identied from their flat power-law, mid-IR SEDs and well- 
characterized optical spectra. Indeed, one of the four source cat- 
egories in our model are type-I quasars and AGNs. Their local 
multi-wavelength LP's and evolution rates are well known, and 
the model predictions have been easily tested on existing data 
(Sect.lOl. 

In contrast, the more numerous type-II AGNs, typically in- 
cluding a mixture of dusty AGN and obscured starburst emis- 
sions, are much more confused among what we called the LIRG 
and ULIRG populations. Because there is no consensus yet on 
the type-II AGN statistical properties, no attempts have been 
made by us to differentiate them from the LIRGs and ULIRGs. 
The broad-band SEDs of our various classes of sources in our 
analysis, reported in Fig.[T7l reveal spectral shapes broadly con- 
sistent with those of star-forming galaxies, with only type-I 
AGNs showing the flat and featureless IR spectra of AGNs. 

From all the above, it is expected that a significant fraction 
of all the IR sources analysed in our paper include some type- 
II AGN contribution. This should be more frequent among our 
ULIRG class and more frequent in mid-IR than far-IR and sub- 
mm samples. For example, Franceschini et al. (2002) suggest 
that assuming approximately 15% of the 24 yum sources to be 
type-II AGN with bright hard X-ray emission would reconcile 
statistical properties of IR and X-ray surveys and the CIRB and 
X-ray cosmological backgrounds. In fact, it is more appropriate 
to think of it in terms of a larger majority of all mid-IR selected 
objects to include AGN components with various degrees of con- 
tribution. 

We can summarize our current understanding of the origin of 
the IR emission in cosmic sources with the statement that star- 
burst and obscured AGN activities are found in many cases to 
occurr simultaneously, with an increased tendency for an AGN 
to appear at increasing luminosity/mass. However, the detailed 
contribution of one or the other energy-production mechanism 
still remains largely undetermined. Our next section will provide 
us with more insight into this problem. 

It should be finally noticed that our previously discussed evo- 
lution pattern of an increased redshift for the peak activity at 
increasing luminosity shown in Fig.l20lcannot in any case be bi- 
ased much by the AGN component, because type-I X-ray AGNs 
have revealed very similar tendencies of their comoving num- 
ber densities vs. redshift as shown here for the whole IR source 
population (Hasinger, Miyaji and Schmidt 2005; Bongiorno et 
al. 2008). 

4.3. Matching the average comoving SFR and the integrated 
stellar mass densities 

Some relevant facts have already emerged from our IR evo- 
lutionary model. The first is that the very fast increase in the 
source emissivity from z = to ~ 1 is happening in a pop- 



ulation of moderate luminosity galaxies (Fig. \T9\ see also Le 
Floc'h et al. 2005), the LIRGs, with typical luminosities of 
Lboi - (1 - 3) X 10" Lq (corresponding to rates of star formation 
of S FR ~ 10-50 M^lyr if dominated by stellar activity, from 
Kennicut's calibration, see below). 

At redshifts higher than ~ 0.7, this population appears to 
stop its evolution and then even reverse it. In parallel, a more 
luminous class of galaxies comes on stage. While providing a 
negligible contribution to the LF's at low-z, this very luminous 
and fast-evolving population overcome the LIRGs at z > 1 . At 
their peak activity epoch around z ~ 1 .5 to 2.5, they show bolo- 
metric luminosities in the ULIRG range {Lboi ^ 10'-Lo,SFRsof 
several hundred M^lyr under the same assumptions as above). 

Altogether, the IR sources manifest a marked bymodal be- 
haviour, particularly evident at z > 1.5. At these redshifts, al- 
though the direct estimates of the LF's available to us from Figs. 
[T4l and [15] do not probe deeply in luminosity and are based 
on moderately numerous samples of a few hundreds sources 
(Rodighiero et al. 2010), the overall shape of the LFs is con- 
strained by our integrated data. Still with reference to the two 
highest-z panels in Fig. [19] on one hand, the pronounced ULIRG 
peak is required to explain the high-z tail in the redshift distri- 
butions of Figs. [T2]and[T3] the high mm number counts in Figs. 
|4]|5] and the observed bright end of the LFs themselves. On the 
other hand, the low-level normalization of the LIRG population 
contribution at lower luminosities and the still flat shape of the 
LFs at these redshifts are constrained by the observed conver- 
gence of the 24 and 70 jim counts and by the requirement not to 
exceed the COBE background. 

At z > 1.5 the bulk of the activity appears to happen in very 
luminous objects: more than 2/3 of it occurs in galaxies more lu- 
minous than 2 lO" Lq. These sites of active SF and obscured ac- 
cretion are presumably also already massive, in consideration of 
the roughly linear proportionality between stellar mass and bolo- 
metric IR emissivity found by various authors (Perez-Gonzales 
et al. 2005; Caputi et al. 2006; Daddi et al. 2007). 

As illustrated in Fig.[T3] our modelling of the statistical prop- 
erties of this luminous population fits those observed for the 
SCUBA galaxies, in terms of both the detailed number counts 
and redshift distributions. Current interpretation identifies these 
ULIRG objects, and the SCUBA sources with the progenitors of 
the massive spheroidal galaxies, where the conversion of large 
amounts of molecular gas into stars are considered to explain 
the huge IR luminosities (e.g. Franceschini et al. 1994; Lilly et 
al. 1999; Blain et al. 2002). Our current modelling of the high-z 
source statistics predicts that only a small fraction of the SCUBA 
objects should be found at z > 3 (see also Rowan-Robinson 
2009). 

Let us now make a working assumption that the bulk of the 
energetics associated with the IR source populations (except of 
course for the few type-I AGNs described in our model) origi- 
nates from young stellar activity. If so, the evolutionary patterns 
in Figs. [18] [19] and l20l would be interpreted to a large extent as 
characters in the evolutionary history of the formation of lumi- 
nous (massive) stars in the universe (massive stars are largely 
responsible for the heating of dust in star-forming regions and 
for dust heating in luminous IR galaxies). We test now this as- 
sumption with the existing complementary information on the 
star-formation history in galaxies based on observations of the 
integrated stellar mass densities. 

To this end, we need to convert our best-fit model distribution 
with redshift of the comoving bolometric IR emissivity in Fig. 
[18] into a comoving rate of stellar formation, in such a way that 
data in this figure can be read in terms of a rate of star formation. 
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Table 1. The stellar mass fraction after gas recycling, as a function of population age for different IMFs of power-law form 
(dn/dlog M oc M""'), and the L/r to SFR calibration 



age 
Gyrs 


Scalo 
(1986) 


Kroupa 
(1993) 


Salpeter 
m = 1.35 


Top-heavy 1 
m = 1.15 


Top-heavy 2 
m = 1.0 


Top-heavy 3 
m = 0.85 


n n 1 
U.Ul 


n no 


U.9/ 


n nn 
U.9U 


n on 


U.o3 


U.oU 


0.05 


0.94 


0.92 


0.85 


0.74 


0.62 


0.53 


0.1 


0.92 


0.90 


0.83 


0.70 


0.58 


0.44 


1.0 


0.86 


0.81 


0.75 


0.61 


0.48 


0.35 


3.0 


0.80 


0.75 


0.72 


0.57 


0.44 


0.32 


7.0 


0.75 


0.72 


0.70 


0.56 


0.42 


0.31 


10. 


0.73 


0.70 


0.65 


0.55 


0.41 


0.30 


15. 


0.71 


0.69 


0.65 


0.54 


0.40 


0.30 


Lis! La 


3.3 10-'" 


2.2 10-'" 


1.63 10-'" 


1.01 10-'" 


0.63 10-'" 


0.46 10-'" 



This conversion is a sensible function of the initial mass func- 
tion (IMF) ruling star-formation. The IMF determines not only 
the relationship between L/r[8 - IQOOfim] and the S FR, but also 
the mass fraction /, of stars remaining at a given galactic epoch, 
considering the fraction of baryonic gas ejected from stars dur- 
ing the evolution and returned to the ISM, or freezed into col- 
lapsed objects, black-holes, white-dwarfs, or neutron stars. The 
fraction /» is obviously a function of the time elapsed from for- 
mation to that of the observation, and depends on the stellar IMF 
(Madau and Pozzetti 2000; Cole et al. 2001). 

To keep our analysis as general as possible, we have made 
use of a set of different IMFs, including significant variations in 
the Salpeter IMF, by increasing or decreasing the relative inci- 
dence of high-mass or low-mass stars. For example, an IMF that 
is more top-heavy than the standard Salpeter one reduces both 
the rate of SF for a given observed Ljr and the average stellar 
residual fraction /» at any galactic epoch. The opposite is ob- 
tained by reducing the relative fraction of high-mass stars (like in 
the Scalo [1986] IMF). To derive the relationship between stel- 
lar mass and bolometric luminosity, we used the Pegase model 
(Fioc and Rocca-Volmerange 1997) by assuming constant SFR 
during 100 Myrs and Category-1 stellar isochrones (Bressan et 
al. 1993). We also assumed that the whole UV flux from young 
stars is absorbed by dust and re-radiated in the IR, a situation 
appropriate for LIRGs and ULIRGs (while it is less clear for 
normal spirals, which do not, however, contribute significantly 
to our energetics). For the standard Salpeter's (1955) power-law 
form for the stellar initial mass function, dn/dlogM oc M^"' for 
stars between 0.1 and 120 Mq, the calibration between the bolo- 
metric luminosity L/r[8 - IQQQfim] and the SFR is 

SFR [MJyr] = K L/r[8 - lOOOjumJ/Lo, (5) 

with 

K =^ 1.63 10"'°, 

consistent with Kennicutt's (1998) reported value. We have also 
considered IMF power-law functions with different slopes m, 
in particular flatter, top-heavier IMF's, for which the constant 
in the relation |5] can change quite significantly. Calibration val- 
ues corresponding to different IMF's are reported in the last row 
of Table [1] together with the stellar residual fractions /»(?) as a 
function of the stellar population age. Data for the Salpeter mass 
function, mass functions of power-law shape between 0.1 and 
120 Mq (dn/dlogM oc M"'") with larger fractions of massive 
stars (Top-heavy 1, 2, 3, corresponding to m - 1.15, 1,0.85), as 
wen as for the Scalo (1986) and the Kroupa et al. (1993) IMFs, 
appear in the table. Going from the Scalo IMF to the top-heavier 



ones, the stellar residual fraction /»(f) and the S FR calibration 
decrease systematically. 

From Eq. |5] a comoving stellar mass density produced by 
these IR star-forming galaxies (type-I AGNs excluded) can be 
easily computed as 

P.,ari> Z) = dz'psFR(z') ■ (^j ■ fMz) " f(z')], (6) 

where Psiar{> z) is the integrated comoving stellar mass func- 
tion at redshift z, ps fr the comoving star-formation rate density, 
scaled from p/R (e.g. as in Fig.fTSI) from eq|5] and where 



dz i/o ■ (1 + z) ■ V(l + zKl + Q„,z) - z(2 + z)Qa) 

is the differential of cosmic time as a function of redshift. The 
function /,[f(z) - t(z')] corrects from the mass converted into 
stars at redshift z' to the present stellar mass at the time t{z). The 
/, values in Table [1] account also for the fraction of baryonic gas 
condensed into collapsed objects. 

4.4. Constraining the history of star formation and AGN 
accretion 

We report in Fig. |2T| predictions by the best-fit IR-evolution 
model concerning the integrated stellar mass akeady assembled 
by the various populations as a function of redshift and compare 
them to the observed integrated stellar mass densities. The three 
predictions correspond to three different adopted IMFs: the stan- 
dard Salpeter (m = 1.35), the Top-heavy 1, and the Kroupa et al. 
(1993) IMFs. Following our analysis in Sect. 14.21 in this calcula- 
tion we have assumed that 80% of the IR emissivity by LIRGs, 
and 70% of that by ULIRGs are due to star formation, with the 
rest of the radiant energy coming from an obscured AGN. This 
corresponds to adopting an AGN fractional contribution to the 
luminosity density of /acn - 0.3 for the ULIRG sources. 

These predictions for the cumulative stellar mass density are 
compared with data derived by Perez-Gonzalez et al. (2008), 
Fontana et al. (2006), and Marchesini et al. (2009) based on a 
large sample of galaxies selected from deep Spitzer/IRAC sur- 
veys from 3.6 to 8 fim. The IRAC data are essential to achieving 
good accuracy in the estimate of the stellar mass content in high- 
redshift galaxies (see also Bundy et al. 2006; Franceschini et al. 
2006; Fontana et al. 2006, among others). We have neglected 
here the (small) correction factors to the observed cumulative 
stellar mass densities in Fig. ISTI corresponding to a change in 
IMF, because the M/L ratio used for the stellar mass estimate 
based on IRAC fluxes is weakly dependent on m. 
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Fig. 22. The stellar mass produced by IR starburst galaxies, according to our reference model, in vaiious redshift bins. Units are 
solar masses per comoving cubic Mpc. Line types as in Fig. [18] Lower panel: the cumulative distribution of stellar mass as a 
function of redshift. The red squared datapoints are estimates of the integrated stellar mass density by Perez-Gonsalez et al. (2008) 
based on deep Spitzer/IRAC surveys, while the cyan data are from Fontana et al. (2006; as indicated by the authors, we applied a 
25% correction to these data to make their mass determination consistent with others in the literature). Blue filled circles are from 
Marchesini et al. (2009). The black line histogram corresponds to our adopted best-fit model for the SF history (in cyan short- 
dashed and red long-dashed those corresponding to our LIRGs and ULIRGs). Upper panel: instantaneous stellar mass produced 
per redshift bin. The normalization factors adopted here are 80% of emitted energy by stars in LIRGs and 70% by ULIRGs (no 
correction for gas recycling by stars). 



As we can see, the three adopted IMFs correspond to sig- 
nificantly different predictions for the integrated stellar mass in 
galaxies as a function of r edshiftB A best fit to the data is clearly 
achieved with the Salpeter IMF, while the Kroupa and Top-heavy 
1 IMFs appear inconsistent with the data, producing systemati- 
cally higher and lower stellar mass distributions as a function of 
redshift. Stronger departures from the Salpeter's function, like 
the Scalo (1986) and the top-heavy IMFs 2 and 3 in Table [1] 
bring into even more serious conflicts with the data. 

There are essentially two sources of uncertainty in our analy- 
sis. The first one is related with the currently still uncertain evo- 
lution of the bolometric luminosity density at z > 2.5, see Sect. 
14.11 a higher emissivity at these redshifts, still partially consis- 
tent with data in our hand, would increase the stellar mass den- 



' Of course, changing the two parameters f^oN and the average stellar 
residual fraction /» does not alter our best-fit modelling of the evolution- 
ary IR emissivity in any way (e.g. in Fig.llSt. but only our interpretation 
of the primary power sources producing it. 



sity at high-z, and consequently also slightly rise the cumulative 
density at lower-z. 

The second uncertainty comes from the poorly known con- 
tribution of AGN gravitational accretion to the observed IR 
bolometric emissivity at high redshifts. The AGN fractions as- 
sumed here are the closest to the data that we discussed in Sect. 
14.21 However, a precise estimate of the contribution of a nu- 
clear non-thermal power source to the bolometric light is still 
not achievable with current data, particularly for high-redshift 
objects (while future studies by Herschel, JWST, and ALMA, 
among others, will off'er further means to disentangle stellar 
and AGN contributions during the evolution of cosmological 
sources). 

In any case, a moderately Top-heavy 1 IMF, like that corre- 
sponding to the dotted red line in Fig. [21] cannot be made en- 
tirely consistent with the data in the figure by decreasing the 
AGN fraction, because our currently assumed value is already 
close to the minimum (it is now < 30%). Perhaps this slightly 
top-heavy IMF could be marginally reconciled with the data if 
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Fig. 21. Predictions of the cumulative distribution of the comov- 
ing stellar mass density as a function of redshift. Units are solar 
masses per comoving cubic Mpc. The black continuous-line his- 
togram corresponds to our adopted model for the SF history and 
a standard Salpeter IMF, the red dotted-line histogram to a Top- 
heavy 1 IMF, the green dashed line to a Kroupa et al. IMF (see 
text for details). The red squared datapoints are estimates of the 
integrated stellar mass density by Perez-Gonzalez et al. (2008) 
based on deep Spitzer/IRAC surveys, while the cyan data are 
from Fontana et al. (2006; as indicated by the authors, we ap- 
plied a 25% correction to these data to make their mass determi- 
nation consistent with others in the literature). Blue filled circles 
are from Marchesini et al. (2009). 



we consider the additional contribution - not accounted for in 
our analysis - of star formation in UV-selected galaxies at high 
redshifts. The contribution of the latter to the bolometric light, 
however, is estimated to be only of the order of 20% of the bolo- 
metric emission (e.g. Hopkins and Beacon 2006; see also Madau 
et al. 1996). 

Also, our prediction based on the Kroupa IMF is difficult 
to reconcile with the data, because it would require adopting a 
very large AGN fraction even for moderate- and low-luminosity 
sources, not justified by the present data. Indeed, the detailed 
shape of the dashed line in Fig. [21] with a larger excess for the 
stellar density at lower-z, would require a higher AGN fraction 
in lower luminosity galaxies, contrary to any expectations. Note, 
finally, that the inconsistency of the dashed-line solution with the 
data would remain large even when adopting a Salpeter IMF for 
the LIRG and using a Kroupa IMF for the ULIRGs. 

Altogether, our comparison of the source bolometric lumi- 
nosity density with the observed stellar mass density as a func- 
tion of redshift reveals that these data are entirely consistent with 
the adoption of a universal Salpeter-like IMF, and contributions 
by AGN accretion to the IR energetics around 20 to 30%. It is 
interesting to stress that this result appears somewhat at vari- 
ance with earlier findings (e.g. Hopkins and Beacon 2006; Perez- 
Gonzalez et al. 2008). These authors have indicated a possible 
redshift- (or luminosity-) dependence of the average stellar IMF 
to reconcile the observational data on the stellar mass and lumi- 
nosity densities. According to them, a top-heavy IMF would be 



required at z > 1 .5 to explain their very high estimated bolomet- 
ric source emissivity at those redshifts. The lower high-z emis- 
sivity indicated by our analysis can be made consistent more 
easily with a universal (Salpeter-like) IMF and standard AGN 
fractions. 

In Fig. |22] we provide further details about our best-guess 
fit to the comoving stellar mass density distribution, with the 
separate LIRG and ULIRG contributions (lower panel), and the 
instantaneous stellar mass produced per redshift bin with no cor- 
rection for gas recycling by stars. As expected, the most lu- 
minous ULIRG population dominates the stellar production at 
z > 1, while the lower luminosity starbursts do it at lower 
redshifts. Non-evolving spirals contribute significantly only at 
z < 0.5 compared to the other populations. 

A substantial AGN fraction among luminous IR MIPS 
sources at high redshifts has been reported by many authors, 
based on hard X-ray observations, population stacking X-ray 
measurements and and IRS spectroscopy (Yan et al. 2006; Sajina 
et al. 2007; Daddi et al. 2007; PoUetta et al. 2008; Martinez- 
Sansigre et al. 2008, Watabe et al. 2009; Fiore et al. 2008, 2009, 
among many others). X-ray stacking by Fiore et al. shows that as 
many as 67% of MIPS objects with ^24 > 550 yuJy and z > 0.6 
contain an AGN, many of which are Compton-thick. Values con- 
sistent with this are inferred by Watabe et al. (2009) from IRS 
mid-IR spectroscopy. These are certainly to be considered as up- 
per limits to the fraction of high-z galaxies containing trace AGN 
emission. In addition, none of these studies can reliably con- 
strain the fractional contribution of the AGN to the bolometric 
source emission. Our completely independent analysis suggests 
that, although a substantial fraction high-z IR sources might con- 
tain an AGN (e.g. Daddi et al. 2007), the AGN contribution 
(from gravitational accretion) to the bolometric flux should be 
small on statistical grounds (< 30%). 

Local remnants of the high-redshift activity revealed by IR 
data would be expected to be massive spheroidal galaxies and 
super-massive black-holes in their centers. The stellar compo- 
nent would make around 30% of the local stellar mass den- 
sity in the universe on average, but would be expected to con- 
centrate on high-mass systems and high-density environment, 
as a consequence of their highly biased formation at high red- 
shifts. Preliminary angular correlation analyses support this view 
(Magliocchetti et al. 2006), but more extensive surveys with 
Herschel will be needed to prove it. 

If we consider, in addition, the close similarity of the evo- 
lutionary patterns as a function of source luminosity described 
in Fig.|20]for the IR-luminous sources with those emerged for 
the type-I X-ray AGNs in Hasinger, Miyaji, and Schmidt (2005) 
and Bongiorno et al. (2008), all this strongly supports the view 
of a tight link between SF and AGN activity in cosmic sources. 
Only high-resolution investigations with future instrumentation 
will reveal whether this is only a statistical connection or the 
signature of a deeper physical relationship among these forming 
structures, such as a ruling effect of one on the other (e.g. AGN 
feedback on the circum nuclear SF). 

As a word of caution, our results from the comparison of 
the stellar mass and SFR densities are based on a combination 
of mid-IR (mostly 24 and 70 jjm) data on the redshift-dependent 
luminosity functions, together with far-IR and sub-mm measure- 
ments of integrated quantities, like number counts and observa- 
tions of the CIRB background radiation (spectral intensity and 
cell-to-cell fluctuations). Although the amount of current data 
strongly constrain our model solutions, we still lack a direct 
evaluation of the bolometric luminosity functions and IR emis- 
sivity of cosmic sources as a function of the cosmic time. The 
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latter will be attempted soon on the basis of the faint far-IR and 
sub-mm source samples that the Herschel Space Observatory is 
starting to produce. Herschel will provide us with a unique op- 
portunity to largely resolve some remaining ambiguities about 
the history of energy production by cosmic sources and its inter- 
pretation. Data from its sensitive photometric imaging cameras, 
PACS and SPIRE, cover a wide wavelength interval between 70 
and 500 ;um, where the bulk of the CIRB energy resides. The 
large systematic programme of cosmological surveys planned 
with Herschel over the whole far-lR and sub-mm spectral region 
(e.g. Griffin et al. 2007; 2010) will soon both resolve the CIRB 
background into its individual sources and, in combination with 
pre-existing optical and Spitzer data over the survey areas, allow 
us to extensively sample the multi-wavelength luminosity func- 
tions. X-ray and radio data will offer additional diagnostic tools 
to disentangle the contributions of the two energy production 
mechanisms. 

5. Conclusions 

We report on our analysis of a large IR database on high- 
redshift galaxies at long wavelengths, the bulk of which we ob- 
tained from recently completed surveys with the Spitzer Space 
Telescope, including unpublished results. Spitzer data are com- 
plemented by sub-mm surveys from the balloon experiment 
BLAST, preliminary results of the Herschel observatory, milli- 
metric observations with large ground-based telescopes, as well 
as diffuse extragalactic background measurements, mostly from 
COBE. 

With these data in hand, we have substatially updated our 
previous (AF2001) modelling of galaxy evolution. Our new re- 
sults agree with AF2001 at z < 1, but much improve their anal- 
ysis at redshifts > 1 , where a very luminous source population 
dominating the IR activity has been established by the new data. 

Our main results are summarized in the following. 

1. The current data completely confirm earlier indications of a 
very rapid increase in galaxy long-wavelength volume emis- 
sivity with redshift up to z ^ 1, paralleled by an increased 
incidence of dust extinction and thermal dust reprocessing in 
high-redshift sources, with respect to locally observed galax- 
ies. This is the fastest evolution rate observed for galaxies 
at any wavelengths [p{z) ^ {I + zY" if averaged over the 
whole galaxy population]. Confusion-limited number counts 
at longer wavelengths (from 70 to 1 100 fim) completely sup- 
port these findings. 

2. All the present data require that the fast evolution observed 
from z-Q to 1 flattens around redshift 1 and keeps approxi- 
mately flat, at least up to z ^ 2.5. At higher redshifts, the co- 
moving emissivity may be required to decrease to avoid ex- 
ceeding the COBE observed intensity of the CIRB, but more 
data from Herschel will be needed for more definite conclu- 
sions. In any case, there is a clear tendency in our redshift- 
dependent bolometric comoving energy density (Fig. [TSl l 
to stay lower than various published results, particularly at 
higher redshifts than z ~ 1.5, based on either large extinc- 
tion corrections, or large spectral extrapolations. 

3. Our exploration of the high-redshift (z > 1) universe in the 
far-IR and sub-mm has provided evidence of a population 
of very luminous galaxies becoming dominant at z > 1. 
The comoving emissivity of the latter is maximal at z - 2 
and shows a much faster convergence with cosmic time than 
lower luminosity systems, whose maximal activity is set 
around z - 1. Then an earlier phase of formation for the 



most luminous/massive galaxies/AGNs is indicated and pro- 
vides supporting evidence of the peculiar evolutionary pat- 
tern named downsizing. At the highest redshifts probed by 
the present analysis (1.5 < z < 3), the luminosity func- 
tions of galaxies appear to be dominated by ultra-luminous 
sources (Ltoi - 10'2 Lq), while the comoving number den- 
sity of lower luminosity galaxies is low and their LFs flat 
(consistent with being not steeper than in the local universe). 
These flat shapes of the LFs at any redshifts indicated by the 
present analysis (Fig.[T9ll might require significant tuning or 
modifications of galaxy formation models. 

4. The IR emissions by cosmic sources analysed here include 
the contributions of both massive-star formation and dust- 
obscured AGN accretion. Many attempts have been made to 
disentangle the two, mostly based on deep hard X-ray imag- 
ing and Spitzer IRS spectroscopy. All these analyses, how- 
ever, remain essentially inconclusive in so far that, in high- 
column-density media, the hard X-ray flux is completely re- 
moved by Compton scattering and the mid-IR hot-dust AGN 
emission may become optically thick. To gain more insight, 
we followed a complementary road of comparing our results 
on the comoving IR emissivity of sources with recent es- 
timates of the redshift-dependent stellar mass functions of 
galaxies based on deep (Spitzer/IRAC) near-IR surveys. In 
our analysis, the comparison of the istantaneous SFR with 
the integrated mass assembled in stars relies on two model 
parameters: the AGN fractional contribution to the bolomet- 
ric IR emission, /agn, and the slope m of the stellar IMF. 
Unlike photonic remnants like the background radiation, suf- 
fering a (1 + z)"' redshift penalty factor, the stellar mass den- 
sity is an unbiased tracer of the past SF activity at the high 
redshifts. 

5. Using our best-fit evolutionary model, we find that the ob- 
served redshift-dependent bolometric source emissivity p/s 
can be naturally reconciled with the observed galaxy mass 
functions by adopting a universal Salpeter IMF (m = 1.35) 
and standard AGN fractions /acw ~ 20 - 30% in LIRGs and 
ULIRGs at all redshifts. Previous analyses referred to higher 
estimated P/r(z) at z > 1.5, hence suggested top-heavy IMFs 
for the more luminous galaxies at higher redshifts or higher 
AGN fractions. 

6. We finally caution that with current data we cannot probe 
deeply into the luminosity functions of high-redshift far-IR 
sources. Statistical constraints were included in our com- 
bined dataset, in particular upper limits on LFs at the faint 
ends from the use of number counts and the CIRB inten- 
sity (including those from cell-to-cell fluctuations). Average 
spectral shapes for populations of sources were derived 
based on the multi-wavelength information from the mid-IR 
to the millimetre. Cosmological surveys with the Herschel 
observatory will soon provide us with data to resolve the 
present residual uncertainties about the history of energy 
production by cosmic sources and its interpretation. 
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Appendix A: Contribution of cosmic sources to tlie 
bacltground radiation: the formalism 

A simple formalism relates background intensity and cell-to-cell 
anisotropics to the statistical properties (luminosity functions 
and number counts) of the contributing sources. 

A. 1 . Source contribution to tlie background intensity 

The differential number counts (sources/unit flux interval/unit 
solid angle) at a given flux S write as 



dN 
dS 



r 



dz — — [L{S ,z),z\ 

dz dS dL 



(A.l) 



where ^[L(5,z),z] is the epoch-dependent (comoving) lumi- 
nosity function (see e.g. Figs. [141 [15] and[T6b and dV/dz is the 
differential (comoving) volume element. Flux S and rest-frame 
luminosity L are related by 



UyK{L,Z) 

Audi ' 



where di is the luminosity distance and K{L,z) - (1 + 
^-jL[v^i+zM jjjg K-correction. The contribution of unresolved 
sources (sources fainter than the detection limit 5 ^[vo] at the fre- 
quency vo) to the background intensity is given by 



>[vo(l+--),.'] 



(A.3) 



(l+Z)X[(l+z)3n,„+nA]l/2' 



for the rather general case of a flat universe with Q„, -i- Qa = 1 . 
Here, j[vo(l + z),z] is the redshift-dependent galaxy comoving 
volume emissivity: 



;[vo(l +z),z] 



inin[Lm„,i.(S<(,z)] 



d\ogLL'^(L,z)K(L,z),(AA) 
dL 



where L^in and L^^x are the minimum and the maximum source 
luminosities. From Eq.(A.3) we can see that, when the counts 
converge like dN/dS oc S^^ or flatter, the contribution by faint 
sources to the background intensity almost becomes insensitive 
to the source minimum flux [/ oc Z«(5„„„) or less]. This property 
has been used by Madau & Pozzetti (2000) to estimate the opti- 
cal background intensity from ultra-deep HST counts of galax- 
ies, by exploiting the convergence of the optical counts fainter 
than niAB ~ 22. Similar properties of convergence of the faint 
IR source counts have been used by Franceschini et al. (2008) to 
estimate the contribution of IR galaxies to the CIRB. 



A.2. Smaii scaie intensity fiuctuations 

In addition to the average integrated flux by all sources in a sky 
area, the background radiation also contains spatial information, 
the cell-to-cell fluctuations, which can be used to further con- 
strain the source flux distribution and spatial correlation proper- 
ties. The usually most important contribution to the cell-to-cell 
intensity fluctuations comes from the stochastic nature of the 
spatial distribution of sources among elementary cells with an 
eff'ective solid angle oj^ff (see below). They can be expressed as 



(A.5) 



What is really measured, however, is not the flux S but the de- 
tector's response x - /(??, (p), f{§, (p) being the angular power 
pattern of the detector Let 



R(x). 



f 

Jo 



dN [x/f(&, <p)] Ids ■ dcjlf{§, ifi) 



(A.6) 



be the mean number of source responses of intensity x in the 
beam. For a Poisson distribution of the number of sources 
producing a response x, its variance equals the mean R(x)dx. 
Adding the variances of all responses up to the cutoff value 
Xc (brighter sources are considered to be individually detected) 
gives the contribution of unresolved sources to fluctuations: 



^ — \ x^ R(x) dx. 
Jo 



(A.7) 



The cutoff Xc is chosen to be some factor q times (51)^; usually 
q = 3-5. The rms background fluctuations (61) imply a sky noise 

9 1/2 

Cconf - {(SI) } for observations with spatial resolution tOe/f- 
(A. 2) Assuming a Gaussian fit with FWHM - 6q to the image's PSF: 



f{6) = exp[-4ie/eofln(2)] 



(A.8) 



with 6o ^ 1.02/i/D to represent the Airy function for diffraction- 
limited observations, we have 



^(xc) — I dxx^R(x) - 
Jo 



(A.9) 



(A. 10) 



— I dxx^ I 
Jo Jo 



dt//N[x/f{t//)]exp(4iffln(2)) 



with iff = 0/6o. In practice, it is custom to compute the confusion 
limit (Tconf from the relation 



00 = cr/ ylnl.{Xc) 



(A.ll) 



by assuming cr - qxc and ^ = 3 or 5 (e.g. Franceschini et al. 
1989). The confusion limit (at 3 or 5 times the rms confusion 
noise (Tconf - c) is then easily computed with the zero's of rela- 
tion lATTI 

The integrated signal D recorded by the detector is the sum 
of the responses x due to all sources in the angular resolution 
element. Its probability distribution function P{D) is informa- 
tive on the amplitude and slope of counts of unresolved sources. 
Scheuer (1957) has shown that, under the assumption of a ran- 
dom sky distribution of sources, its Fourier transform, pito), is a 
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PSF FWHM 00 (arcsec) 



Fig. A.l. Estimated 3cr confusion limits as a function of the FWHM angle of a Gaussian fit to the PSF (see Appendix A and B), 
0' = 3, for observations at various effective wavelengths. The latter are reported with the following colour and line codes: black dots: 
A =35 jum; red continuous: 70 fim; cyan short dash: 100 //m; green dots: 160 fim; black continuous: 250 fim; red long dash: 500 fim; 
black short dash: 850 yum. 



simple function of the FT r(oj) of R(x): piu) - exp[r((y) - r(0)]. 
It follows: 



p{u) &xp{-2m(jjD) doj - 

> 

exp [r((jj) - r(0) - InicoD] dco - 
2 ^ exp |- ^ R(x) [1 - cos(27rwx)] c/xj ■ 

■ cos I R{x) sin(2;rwx) dx - IttojD 
Jo 



(A. 12) 
(A.13) 
(A. 14) 



da>. 



This synthetic P(D) has to be convolved with the (typically 
Gaussian) distribution of the instrumental noise to be compared 
with the observations. Rather than in terms of the total flux per 
sky beam D, it is more frequent to refer the fluctuation analysis 
to the flux per unit sky area Da - D/a>eff, where a>eff is the 
imager effective area and is defined as 



J^' xR(x)dx 
/(< X,) 



(A. 15) 



with /(< Xc) being the residual background, from Eq. iA.3\ . It 
turns out that 



/(??)£/□ = ;r J d9\f{0) =7T{0o/2f/ln(2) (A.16) 



for our Gaussian representation of the PSF in eq. (IA.8b . 

Assumed that the number count distribution below the de- 
tection limit can be represented as a power-law, A^(> S) = 
K(S ISkT^^, then eq. IIA.14II can be integrated to get (Condon 
1974): 



^conf ~ 



q 



2-/5 



J^/()?, ipfd^l 



(A. 17) 
(A.18) 



which allows estimation of the slope of the counts (fi) below 
the detection limit from a given measured value of the cell-to- 
cell fluctuations cTconf- This constraint on NiS) applies down to 
a flux limit corresponding to ~ 1 source/beam. Assumed that 
S k represents the confusion limit {Sk = q x crc„nf) of a survey 
having an areal resolution oJeff, then eq. lA.18l further simplifies 
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to a relation between the number of sources K resolved by the 
survey (and brighter than S k) and the parameters q and /5: 



Table B.l. Euclidean-normalized differential and integral num- 
ber counts at A^ff = 35 and 70 fj.m by our best-fit model 



K : 



2-/3 1 



(A.19) 



this implies the confusion limit to occur at the flux correspond- 
ing to an areal density of (/3q^/[2 -jS]) ' sources per unit beam 
area atg-. For Euclidean counts and q = 3, this corresponds to 1 
source/27 beams. 



Appendix B: Model predictions for future 
observations 

B.1. Source number counts 

Predictions for (differential and integral) source number counts 
are reported in Tables lBn to lB.4l based on our best-fit model. 

In each table the left column reports the logarithm of 
flux, and the second and third columns give the (Euclidean- 
normalized) differential and integral total number counts for 
the first effective wavelength indicated in the table caption. The 
fourth and fifth columns give the differential and integral counts 
for the second wavelength in the caption. 

The units for the differential counts are the number of 
sources per square degree per unit flux interval in mJy, divided 
by the Euclidean function S^^^^. Those for the integral counts 
are number per square degree. 

B.2. Predictions for source confusion limitations 

We calculated predictions for the source confusion limits based 
on the best-fit model for observations at various wavelengths us- 
ing relation lA. 1 ll in Appendix A and calculating the zero's of the 
relation as a function of the FWHM Go of the observation at a 
given wavelength. This definition of confusion is sometimes re- 
ferred to as photometric confusion, dealing with the photometric 
noise induced by source pile-up in the map's elementary beam. 

As mentioned in Appendix A, for a diffraction-limited imag- 
ing survey the FWHM angle 6q of a Gaussian fit to the PSF Airy 
function is given by 



00 



1.02^ 
D 



(B.l) 



with A being the effective observational wavelength and D the 
primary collector diameter If the imager is not diffraction- 
limited, 0() is in any case the PSF FWHM. We report in Fig. lA.ll 
our estimated 3cr confusion limits for various far-IR and sub-mm 
wavelengths assuming a value q -!>. 
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3.046 


-0.2809 


4.179 


0.8207 


2.3 


3.047 


-0.5711 


4.154 


0.4974 


2.5 


3.051 


-0.8624 


4.131 


0.1821 


2.7 


3.052 


-1.176 


4.114 


-0.1331 


2.9 


3.057 


-1.481 


4.099 


-0.4463 


3.1 


3.052 


-1.781 


4.086 


-0.7606 


3.3 


3.045 


-2.091 


4.072 


-1.079 


3.5 


3.042 


-2.406 


4.056 


-1.401 


3.7 


3.032 


-2.72 


4.034 


-1.73 


3.9 


3.01 


-3.056 


4.006 


-2.068 
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Table B.2. Euclidean-normalized differential and integral num- Table B.3. Euclidean-normalized differential and integral num- 
ber counts at Aef/ = 100 and 160 //m by our best-fit model ber counts at Ae/f = 250 and 350 /zm by our best-fit model 



log(S) 


log[dN/dS] 


logm> S)] 


log[dN/dS] 


logm>S)} 


log(S) 


log[dNldS] 


log[N{>S)-\ 


logidNldS} 


log[N{>^S)] 


mJy 


deg-'mjy-^ 


deg-' 


deg~'mJy'^ 


deg-' 


mJy 


deg-'mjy-^ 


deg-' 


deg-'mjy-^ 


deg-' 


-1.9 


1.811 


5.171 


1.872 


5.307 


-1.9 


1.866 


5.392 


1.835 


5.404 


-1.7 


2.074 


5.111 


2.154 


5.256 


-1.7 


2.188 


5.348 


2.188 


5.364 


-1.5 


2.335 


5.047 


2.429 


5.202 


-1.5 


2.493 


5.298 


2.522 


5.313 


-1.3 


2.593 


4.981 


2.698 


5.144 


-1.3 


2.785 


5.242 


2.834 


5.253 


-1.1 


2.853 


4.911 


2.965 


5.081 


-1.1 


3.065 


5.181 


3.126 


5.182 


-0.9 


3.115 


4.836 


3.23 


5.015 


-0.9 


3.339 


5.113 


3.401 


5.102 


-0.7 


3.375 


4.752 


3.494 


4.945 


-0.7 


3.61 


5.038 


3.663 


5.013 


-0.5 


3.633 


4.658 


3.761 


4.867 


-0.5 


3.876 


4.956 


3.912 


4.91 


-0.3 


3.882 


4.549 


4.027 


4.78 


-0.3 


4.138 


4.862 


4.147 


4.794 


-0.1 


4.118 


4.423 


4.287 


4.678 


-0.1 


4.392 


4.754 


4.369 


4.663 


0.1 


4.332 


4.275 


4.536 


4.558 


0.1 


4.631 


4.624 


4.576 


4.514 


0.3 


4.52 


4.099 


4.765 


4.415 


0.3 


4.847 


4.469 


4.761 


4.344 


0.5 


4.674 


3.891 


4.963 


4.24 


0.5 


5.032 


4.282 


4.917 


4.139 


0.7 


4.788 


3.649 


5.121 


4.025 


0.7 


5.172 


4.054 


5.04 


3.896 


0.9 


4.859 


3.372 


5.229 


3.762 


0.9 


5.261 


3.774 


5.118 


3.614 


1.1 


4.887 


3.055 


5.277 


3.45 


1.1 


5.297 


3.442 


5.139 


3.279 


1.3 


4.881 


2.714 


5.266 


3.084 


1.3 


5.267 


3.065 


5.085 


2.873 


1.5 


4.847 


2.363 


5.207 


2.682 


1.5 


5.166 


2.625 


4.931 


2.363 


1.7 


4.798 


2.007 


5.112 


2.263 


1.7 


5.009 


2.112 


4.696 


1.742 


1.9 


4.744 


1.651 


5.002 


1.848 


1.9 


4.822 


1.575 


4.418 


1.093 


2.1 


4.689 


1.297 


4.891 


1.445 


2.1 


4.626 


1.1 


4.144 


0.5516 


2.3 


4.636 


0.9495 


4.786 


1.051 


2.3 


4.453 


0.6731 


3.92 


0.1081 


2.5 


4.587 


0.6062 


4.692 


0.6677 


2.5 


4.322 


0.2747 


3.769 


-0.2791 


2.7 


4.544 


0.2677 


4.61 


0.2964 


2.7 


4.224 


-0.09666 


3.672 


-0.6411 


2.9 


4.504 


-0.06816 


4.541 


-0.05632 


2.9 


4.149 


-0.4531 


3.6 


-0.9954 


3.1 


4.466 


-0.4033 


4.484 


-0.3998 


3.1 


4.088 


-0.8021 


3.54 


-1.356 


3.3 


4.431 


-0.7385 


4.439 


-0.7369 


3.3 


4.034 


-1.151 


3.48 


-1.716 


3.5 


4.397 


-1.069 


4.398 


-1.074 


3.5 


3.984 


-1.503 


3.424 


-2.071 


3.7 


4.367 


-1.397 


4.362 


-1.412 


3.7 


3.937 


-1.849 


3.375 


-2.415 


3.9 


4.342 


-1.72 


4.327 


-1.748 


3.9 


3.895 


-2.188 


3.33 


-2.758 
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Table B.4. Euclidean-normalized differential and integral num- 
ber counts at Aef/ = 500 and 850 //m by our best-fit model 



log(S) log[dN/dS] log[N(>S)] logWN/dS] log[N(> S)] 



mjy 


deg 'mJy ' 


deg-' 


deg 'mJy ' 


deg-' 


-1.9 


1.859 


5.382 


2.049 


5.214 


-1.7 


2.225 


5.336 


2.346 


5.101 


-1.5 


2.567 


5.276 


2.585 


4.955 


-1.3 


2.882 


5.2 


2.765 


4.775 


-1.1 


3.17 


5.108 


2.877 


4.584 


-0.9 


3.431 


4.998 


2.932 


4.396 


-0.7 


3.66 


4.868 


3.006 


4.251 


-0.5 


3.855 


4.717 


3.157 


A.\6A 


-0.3 


4.017 


4.552 


3.376 


4.097 


-0.1 


4.164 


4.384 


3.644 


3.993 


0.1 


4.322 


4.223 


3.916 


3.808 


0.3 


4.484 


4.071 


4.105 


3.552 


0.5 


4.642 


3.88 


4.161 


3.219 


0.7 


4.781 


3.629 


4.064 


2.759 


0.9 


4.846 


3.312 


3.812 


2.116 


1.1 


4.796 


2.899 


3.425 


1.271 


1.3 


4.607 


2.338 


2.952 


0.4295 


1.5 


4.301 


1.607 


2.491 


-0.3523 


1.7 


3.926 


0.8353 


2.127 


-0.9304 


1.9 


3.558 


0.1836 


1.862 


-1.363 


2.1 


3.266 


-0.3064 


1.686 


-1.769 


2.3 


3.076 


-0.6918 


1.558 


-2.198 


2.5 


2.961 


-1.048 


1.443 


-2.623 


2.7 


2.888 


-1.411 


1.328 


-3.04 


2.9 


2.818 


-1.788 


1.207 


-3.447 


3.1 


2.752 


-2.159 


1.075 


-3.892 


3.3 


2.688 


-2.522 


0.9155 


-4.384 


3.5 


2.629 


-2.87 


0.7131 


-4.916 


3.7 


2.567 


-3.242 


0.4657 


-5.485 


3.9 


2.491 


-3.641 


0.1755 


-6.118 



